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HIGH THROUGHPUT DIRECTED EVOLUTION BY RATIONAL 

MUTAGENESIS 

RELATED APPLICATIONS 

Benefit of priority under 35 U.S.C. § 1 1 9(e) is claimed to U.S. 
provisional application Serial No. 60/31 5,382, filed August 27, 2001, to 
Manuel Vega, Lila Drittanti and Marjorie Flaux, entitled "HIGH THROUGH- 
PUT DIRECTED EVOLUTION BY RATIONAL MUTAGENESIS." The subject 
matter of this application is incorporated in its entirety by reference thereto. 
FIELD OF INVENTION 

Processes and systems for the high throughput directed evolution 
of peptides and proteins, particularly those that act in complex biological 
settings, are provided. The proteins and peptides include, but are not 
limited to, intracellular proteins, messenger/signaling/hormone proteins 
and viral proteins. 
BACKGROUND 

Directed evolution refers to biotechnological processes for 
optimizing the activity of proteins by means of random changes 
introduced into selected respective genes. Directed evolution involves the 
creation of a library of mutated genes, and then selection of the mutants 
that encode proteins having desired properties. The process can be an 
iterative one in which gene products that have improvement in a desired 
property are subjected to further cycles of mutation and screening. 
Directed evolution provides a way to adapt natural proteins to work in 
new chemical or biological environments, and/or to elicit new functions. 
The potential plasticity of proteins is such that chances exist that for 
every new challenge, such as a new environment and desired new or 
altered activity, it should be possible, given a sufficient pool of modified 
proteins (or encoding nucleic acids), that an appropriately 'evolved' 
protein could be found that would have a desired activity. The problem is 
in generating and then identifying the appropriate sequence. 

There have been practical approaches to this problem (see, e.g., 
U.S. Patent Nos. 6,096,548; 6,117,679; 6,165,793; 6,180,406; 
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6,132,970; 6,171,820; 6,238,884; 6,174,673; 6,057,103; 6,001,574; 
5,763,239; 5,837,500; 5,571,698; 6,156,509; 5,723,323; 5,862,514; 
5,871,974; 5,779,434 and others). Typically theses approaches are of 
two types. One is a purely "rational" approach that is based on the 
assumption that the optimized proteins can be rationally designed. This, 
however, requires sufficient information regarding the laws that govern 
protein folding, molecular interactions, intra-molecular forces and other 
dynamics of protein activity. This rational approach is extremely 
dependent on a number of variables and parameters that are not known. 
Consequently, although useful in some specific cases and applications, 
the rational approach intended to 'predict' protein structure remains 
limited in applicability. 

In contrast to the rational approach, random approaches have also 
been employed. One random approach requires synthesis of all possible 
protein sequences or a statistically sufficient large number of proteins and 
then screening them to identify proteins having the desired activity or 
property. Since the resources to synthesize all possible theoretical 
sequences of a single protein is not possible, this approach is 
impracticable. Other random approaches are based on gene shuffling 
methods, which are PCR-based methods that generate random 
rearrangements between two or more sequence-related genes to 
randomly generate variants of the gene. 

The development and scope of directed evolution, thus, has been 
limited, and its potential remains to be exploited. In order to exploit the 
potential of directed evolution, alternative approaches for generating and 
identifying evolved proteins are needed. It is an object herein to provide 
methods and products to exploit the potential of directed evolution. 
SUMMARY 

Provided herein are methods for performing directed evolution for 
the optimization of proteins that function in complex biological settings. 
Methods of high throughput directed evolution of proteins are provided. 
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In practicing the methods, each molecule is individually designed, 
produced, processed, screened and tested in a high throughput format. 
Neither random or combinatorial methods nor mixtures of molecules are 
used. 

The methods provided herein include the steps of identifying a 
protein target of interest; obtaining nucleic acids that encode the target, 
which may be from any source, such as a natural library, a collection 
generated by known gene shuffling techniques and related methods, and, 
then creating variants of the proteins using methods for rational 
mutagenesis provided herein. Whatever method is used to select or 
generate the nucleic acids encoding the protein targets, each molecule is 
processed and screened separately in a high throughput format. 

The nucleic acids encoding each variant are individually screened. 
They can be screened in any suitable assay, including cell-based assays 
and biochemical assays. For cell based assays, each nucleic acid 
molecule is introduced into an expression vector for expression in a 
bacterial cell or into a vector for expression in a eukaryotic host cell. In all 
instances, the nucleic acids of interest are introduced into host cells in an 
expression vector, such as by transfection for bacterial hosts and 
transduction with viral vectors into eukaryotic hosts with viral vectors. 

Each variant is introduced into a host and the resulting cells are 
maintained separately, such as in an addressable array of wells in a 
microtiter plate or other substrate with discrete locations for performing 
reactions or retaining molecules of interest. Typical formats are 96 loci, 
and multiples thereof (384, 1536, 3072, . . . 96 x n, where n is 1 to any 
number desired, such as 10, 20, 30, 50 . . .100), although any 
convenient number of loci may be employed. 

Since the process is conducted in a high throughput format, for 
many embodiments, it is often important to assess the relative numbers 
of transformed, transduced or transfected cells. Hence the relative (or 
actual) titer of the vector, such as the recombinant viral vector, must be 



37851-911 



known to permit analysis of results. For high throughput formats, it is 
important to assess the relative or actual concentration of the viral vector 
(or plasmid) so that results can be compared among all cells and variants. 
Methods for titering (determining the concentration) of the nucleic acid 
encoding the variant and/or the recombinant virus are also provided. 

The processes require accurate titering of the viruses in a collection 
or among collections (libraries) so that the activities of the screened 
mutant proteins can be compared. Provided are general methods for the 
quantitative assessment of the parameters of activity corresponding to 
the individual variants in the library, based upon intracellular serial dilution 
generated by precise tittering with the gene transfer viral vectors Any 
method permits accurate titering may be used, including that described 
nlnternational PCT application No. PCT/FR01/01366, based on French 
application n° 0005852, filed 9 May 2000, and published as International 
PCT application No. WO 01/186291 . A method of titering, designated 
Tagged Replication and Expression Enhancement Technology (TREE") is 
provided herein. 

Each of the different cells is separately screened by a suitable 
assay, and the results analyzed. Methods for assessing the interactions 
in biological systems, such as a Hill-based analysis (see, published 
International PCT application No. WO 01/44809 based on PCT n° 
PCT/FR00/03503, Dec, 2000, and the description herein), or a second 
order polynomial or other algorithm that describes the interaction between 
cells and biological agents to select variants that have a desired property 
are employed in the processes herein. 

A semi-rational method for evolution of proteins that is particularly 
designed for use in the methods herein or in any method that uses 
"evolved" proteins is also provided. The method, which is based on an 
amino-acid scanning protocol, is for rationally designing the variants for 
use in the directed evolution and selection method, and can employ 
iterative processing of the steps of the high throughput methods provided 
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herein. In this method, once the target protein or domain is identified, 
nucleic acid molecules encoding variants are prepared. Each variant 
encoded by the nucleic acid molecules has a single amino acid replaced 
with another selected amino acid, such as alanine (Ala), glycine (Gly), 
serine (Ser) or any other suitable amino acid, typically one selected to 
have a neutral effect on secondary and tertiary structure. The resulting 
series of variants are separately screened in the high throughput format 
provided herein, and those that have a change in the target activity are 
selected and the modified amino acids are designated "hits." Nucleic acid 
molecules encoding proteins in which each hit position is replaced by the 
eighteen remaining amino acids then are synthesized and the resulting 
collection of molecules are screened, such by introduction into host cells, 
and the proteins that result in a improvement of a targeted activity, are 
identified. Such proteins are designated "leads." Leads may be further 
modified by producing proteins that have combinations of the mutations 
identified in the leads. This method, which does not require any 
knowledge of the structure of a target protein, permits precise control of 
locations where changes are introduced and also the amount of change 
that is introduced. 

The high throughput directed evolution processes provided herein 
include the use of virus libraries containing mutant versions of a gene; 
viral libraries of such mutant genes are also provided. 

Reporter cells are infected with the titered viruses that encode the 
mutant genes. The mutant genes are expressed and read-out data from 
either biochemical or cell-based assays, while isolating each mutant/virus 
physically from the others (i.e. one-by-one analysis), are collected and 
analyzed. Serial dilution assays (i.e. a series of dilutions for each 
individual mutant/virus in the library) are used and the biochemical/cell- 
based assays are performed on each single dilution for each individual 
mutant/virus. Analysis of the serial dilution readout-data can be 
performed using any method of analysis that permits one-by-one 



• 



37851-911 



comparisons. Hill-based analysis (see, published International PCT 
application No. WO 01/44809 based on PCT n° PCT/FR00/03503, Dec, 
2000, and the description herein) are employed for analysis of the data. 
Protein/protein domain variants identified using the methods are 
5 also provided. Also provided are nucleic acid molecules and proteins and 
polypeptides produced by the methods and viruses and cells that contain 
the nucleic acid moleucles and proteins. 

In an exemplary embodiment of methods provided herein, the 
process of rational directed evolution provided herein is applied to the 

10 AAV rep gene. The resulting recombinant rep protein variants and rAAV 
are also provided. Among the rep proteins are those that result in 
increased rAAV production in rAAV that encode such mutants, thereby, 
among a variety of advantages, offering a solution to the need in the gene 
therapy industry to increase the production therapeutic vectors without 

15 up-scaling manufacturing. 

Thus, for exemplification, some methods provided herein have been 
used to identify amino acid "hit" positions in adeno-associated virus 
(AAV) rep proteins that are relevant for AAV or rAAV production. Those 
amino acid positions are such that a change in the amino acid leads to a 

20 change in protein activity either to lower activity or to higher activity 
compared to native-sequence Rep proteins. The hit positions were then 
used to generate further mutants designated "leads." Provided herein are 
the resulting mutant rep proteins that result in either higher or lower 
levels of AAV or rAAV virus compared to the wild-type (native) Rep 

25 protein(s). 

In addition to enhancing AAV production, among the rep mutants 
are those that inhibit papillomavirus (PV) and PV-associated diseases, 
including certain cancers and human immunodeficiency virus (HIV) and 
HIV-associated diseases. 
30 Systems and computer controlled systems for performing the high 

throughput processes are also provided. 

-6- 
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DESCRIPTION OF THE FIGURES 

FIGURES 1 summarize various exemplary embodiments of the high 
throughput processes provided herein. FIGURE 1A depicts an 
embodiment of the process in which an amino acid scan is employed to 
generate a library of mutants, which are then introduced into viral 
vectors, such as an adeno-associated viral vector (AAV), an herpes virus, 
such as herpes simplex virus (HSV) and other herpes virus vectors, a 
vaccinia virus vector, retroviral vectors, such as MuMLV, MoMLV, feline 
leukemia virus, and HIV and other lentiviruses, adenovirus vectors and 
other suitable viral vector, each member of the library is individually 
tested and phenotypically characterized to identify HITS. FIGURE 1B 
summarizes round 2 in which LEADS are developed by mutagenesis at 
and/or surrounding the positions identified as HITS; FIGURE 1C 
summarizes the optional next round in which recombination among 
LEADS is performed to further optimize the LEADS; FIGURE 1D depicts 
the process in mammalian cells; and FIGURE 1E depicts the process in 
bacterial cells. 

FIGURE 2A depicts an exemplary titering process (in this instance 
the TREE™ for titering AAV) in a 96 well format; FIGURE 2B shows the 
results and analysis of a titering process performed using the TREE™ 
procedure; and FIGURE 2C shows an exemplary calibration curve for the 
calculation of the titer using the TREE" method. 

FIGURE 3A and 3B depict "HITS" and "LEADS" respectively for 
identification of AAV rep mutants "evolved" for increased activity. 

FIGURE 4 shows the genetic map of AAV, including the location of 
promoters, and transcripts; amino acid 1 of the Rep 78 gene is at 
nucleotide 321 in the AAV-2 genome. 

FIGURES 5A and 5B show the alignment of amino acid sequences 
of Rep78 among AAV-1; AAV-6; AAV-3; AAV-3B; AAV-4; AAV-2; AAV- 
5 sequences, respectively; the hit positions with 100 percent homology 
among the serotypes are bolded italics, where the position is different 



(compared to AAV-2, no. 6 in the Figure) in a particular serotype, it is in 
bold; a sequence indicating relative conservation of sequences among 
the serotypes is labeled "C". 

Legend: 

1 is AAV-1 ; 2 is AAV-6, 3 is AAV-3, 4 is AAV-3B, 
5 is AAV-4, 6 is AAV-2, and 7 is AAV-5; 
"." where the amino acid is present > 20%; 

:" where the amino acid is present > 40%; 
" + " where the amino acid is present > 60%; 
y 10 "*" where the amino acid is present > 80%; and 

Pi where the amino acid is the same amongst all 

p serotypes depicted it is represented by its single letter 

code. 
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DETAILED DESCRIPTION 

15 A. Definitions 

Unless defined otherwise, all technical and scientific terms used 
herein have the same meaning as is commonly understood by one of skill 
in the art to which this invention belongs. All patents, patent 
applications, published applications and publications, Genbank sequences, 

20 websites and other published materials referred to throughout the entire 
disclosure herein are, unless noted otherwise, incorporated by reference 
in their entirety. In the event that there are a plurality of definitions for 
terms herein, those in this section prevail. 

As used herein, directed evolution refers to mehods that adapt" 

25 natural proteins or protein domains to work in new chemical or biological 
environments and/or to elicit new functions. It is more a more broad- 
based technology than DNA shuffling. 

As used herein, high-throughput screening (HTS) refers to 
processes that test a large number of samples, such as samples of test 

30 proteins or cells containing nucleic acids encoding the proteins of interest 
to identify structures of interest or the identify test compounds that 
interact with the variant proteins or cells containing them. HTS 
operations are amenable to automation and are typically computerized to 
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handle sample preparation, assay procedures and the subsequent 
processing of large volumes of data. 

As used herein, DNA shuffling is a PCR-based technology that 
produces random rearrangements between two or more sequence-related 
genes to generate related, although different, variants of given gene. 

As used herein, "hits" are mutant proteins that have an alteration in 
any attribute, chemical, physical or biological property in which such 
alteration is sought. In the methods herein, hits are generally generated 
by systematically replacing each amino acid in a protein or a domain 
thereof with a selected amino acid, typically Alanine, Glycine, Serine or 
any amino acid, as long as each residue is replaced with the same 
residue. Hits may be generated by other methods known to those of skill 
in the art and tested by the high throughput methods herein. For 
purposes herein a Hit typically has activity with respect to the function of 
interest that differs by at least 10%, 20%, 30% or more from the wild 
type or native protein. The desired alteration, which is generally a 
reduction in activity, will depend upon the function or property of interest. 

As used herein, "leads" are "hits" whose activity has been 
optimized for the particular attribute, chemical, physical or biological 
property. In the methods herein, leads are generally produced by 
systematically replacing the hit loci with all remaining 18 amino acids, and 
identifying those among the resulting proteins that have a desired activity. 
The leads may be further optimized by replacement of a plurality of "hit" 
residues. Leads may be generated by other methods known to those of 
skill in the and tested by the high throughput methods herein. For 
purposes herein a lead typically has activity with respect to the function 
of interest that differs from the native activity, by a desired amount and is 
at by at least 10%, 20%, 30% or more from the wild type or native 
protein. Generally a Lead will have an activity that is 2 to 10 or more 
times the native protein for the activity of interest. As with hits, the 



change in the activity is dependent upon the activity that is "evolved." 
The desired alteration will depend upon the function or property of 
interest. 

As used herein, MOI is multiplicity of infection. 

As used herein, ip, with reference to a virus or recombinant vector, 
refers to a titer of infectious particles. 

As used herein, pp refers to the total number of vector (or virus) 
physical particles 

As used herein, biological and pharmacological activity includes any 
activity of a biological pharmaceutical agent and includes, but is not 
limited to, biological efficiency, transduction efficiency, gene/transgene 
expression, differential gene expression and induction activity, titer, 
progeny productivity, toxicity, cytotoxicity, immunogenicity, cell 
proliferation and/or differentiation activity, anti-viral activity, 
morphogenetic activity, teratogenetic activity, pathogenetic activity, 
therapeutic activity, tumor suppressor activity, ontogenetic activity, 
oncogenetic activity, enzymatic activity, pharmacological activity, 
cell/tissue tropism and delivery. 

As used herein, "output signal" refers to parameters that can be 
followed over time and, if desired, quantified. For example, when a virus 
infects a cell, the infected cell undergoes a number of changes. Any 
such change that can be monitored and used to assess infection, is an 
"output signal," and the cell is referred to as a "reporter cell." Output 
signals include, but are not limited to, enzyme activity, fluorescence, 
luminescence, amount of product produced and other such signals. 
Output signals include expression of a viral gene or viral gene product, 
including heterologous genes (transgenes) inserted into the virus. Such 
expression is a function of time ("t") after infection, which in turn is 
related to the amount of virus used to infect the cell, and, hence, the 
concentration of virus ("s") in the infecting composition. For higher 
concentrations the output signal is higher. For any particular 
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concentration, the output signal increases as a function of time until a 
plateau is reached. Output signals may also measure the interaction 
between cells, expressing heterologous genes, and biological agents 

As used herein, adeno-associated virus (AAV) is a defective and 
non-pathogenic parvovirus that requires co-infection with either 
adenovirus or herpes virus for its growth and multiplication, able of 
providing helper functions. A variety of serotypes are known, and 
contemplated herein. Such serotypes include, but are not limited to: 
AAV-1 (Genbank accession no. NC002077; accession no. VR-645); AAV- 
2 (Genbank accession no. NC001401; accession no. VR-680); AAV-3 
(Genbank accession no. NC001729; acession no. VR-681); AAV-3b 
(Genbank accession no. NC001863); AAV-4 (Genbank accession no. 
NC001829; ATCC accession no. VR-646 ); AAV-6 (Genbank accession 
no.NC001729); and avian associated adeno-virus (ATCC accession no. 
VR-1449). The preparation and use of AAVs as vectors for gene 
expression in vitro and for in vivo use for gene therapy is well known 
(see, e.g., U.S. Patent Nos. 4,797,368, 5,139,941, 5,798,390 and 
6,127,175; Tessier eta/. (2001) J. Virol. 75:375-383; Salvetti eta/. 
(1998) Hum Gene Ther 20:695-706; Chadeuf eta/. (2000) J Gene Med 
2:260-268). 

As used herein, the activity of a Rep protein or of a capsid protein 
refers to any biological activity that can be assessed. In particular, 
herein, the activity assessed for the rep proteins is the amount (i.e., titer) 
of AAV produced by a cell. 

As used herein, the Hill equation is a mathematical model that 
relates the concentration of a drug (i.e., test compound or substance) to 
the response being measured 

Waxier: 

y = [D] n + [D 50 ] n 
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where y is the variable being measured, such as a response, signal, y max is 
the maximal response achievable, [D] is the molar concentration of a 
drug, [D 50 ] is the concentration that produces a 50% maximal response to 
the drug, n is the slope parameter, which is 1 if the drug binds to a single 
site and with no cooperativity between or among sites. A Hill plot is log 10 
of the ratio of ligand-occupied receptor to free receptor vs. log [D] (M). 
The slope is n, where a slope of greater than 1 indicates cooperativity 
among binding sites, and a slope of less than 1 can indicate heterogeneity 
of binding. This general equation has been employed for assessing 
interactions in complex biological systems (see, published International 
PCT application No. WO 01/44809 based on PCT n° PCT/FR00/03503, 
see, also, EXAMPLES). 

As used herein, in the Hill-based analysis (published International 
PCT application No. WO 01/44809 based on PCT n° PCT/FR00/03503), 
the parameters, rr,K,T,€,fi,0, are as follows: 

n potency of the biological agent acting on the assay (cell- 
based) system; 

k constant of resistance of the assay system to elicit a 
response to a biological agent; 

e is global efficiency of the process or reaction triggered by the 
biological agent on the assay system; 

r is the apparent titer of the biological agent; 
6 is the absolute titer of the biological agent; and 
q is the heterogeneity of the biological process or reaction. 
In particular, as used herein, the parameters n (potency) or k 
(constant of resistance) are used to respectively assess the potency of a 
test agent to produce a response in an assay system and the resistance 
of the assay system to respond to the agent. 

As used herein, ^(efficiency), is the slope at the inflexion point of 
the Hill curve (or, in general, of any other sigmoidal or linear 
approximation), to asses the efficiency of the global reaction (the 
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biological agent and the assay system taken together) to elicit the 
biological or pharmacological response. 

As used herein, r (apparent titer) is used to measure the limiting 
dilution or the apparent titer of the biological agent. 

As used herein, 0 (absolute titer), is used to measure the absolute 
limiting dilution or titer of the biological agent. 

As used herein, q (heterogeneity) measures the existence of 
discontinuous phases along the global reaction, which is reflected by an 
abrupt change in the value of the Hill coefficient or in the constant of 
resistance. 

As used herein, a library of mutants refers to a collection of 
plasmids or other vehicles that carrying (encoding) the gene variants, 
such that individual plasmid or other vehicles carry individual gene 
variants. When a library of proteins is contemplated, it will be so-stated. 

As used herein, a "reporter cell" is the cell that "reports", i.e., 
undergoes the change, in response to introduction of the nucleic acid 
infection and, therefore, it is named here a reporter cell. 

As used herein, "reporter" or "reporter moiety" refers to any moiety 
that allows for the detection of a molecule of interest, such as a protein 
expressed by a cell. Typical reporter moieties include, for example, 
fluorescent proteins, such as red, blue and green fluorescent proteins. 
For expression in cells, nucleic acid encoding the reporter moiety can be 
expressed as a fusion protein with a protein of interest or under to the 
control of a promoter of interest. 

As used herein, a titering virus increases or decreases the output 
signal from a reporter virus, which is a virus that can be detected, such 
as by a detectable label or signal. 

As used herein, phenotype refers to the physical or other 
manifestation of a genotype (a sequence of a gene). In the methods 
herein, phenotypes that result from alteration of a genotype are assessed. 
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As used herein, activity refers to the function or property to be 
evolved An active site refers to a site(s) responsible or that participates 
in conferring the activity or function. The activity or active site evolved 
(the function or property and the site conferring or participating in 
conferring the activity) may have nothing to do with natural activities of 
a protein. For example, it could be an 'active site' for conferring 
immunogenicity (immunogenic sites or epitopes) on a protein. 

As used herein, the amino acids, which occur in the various amino 
acid sequences appearing herein, are identified according to their known, 
three-letter or one-letter abbreviations (see, Table 1). The nucleotides, 
which occur in the various nucleic acid fragments, are designated with 
the standard single-letter designations used routinely in the art. 

As used herein, amino acid residue refers to an amino acid formed 
upon chemical digestion (hydrolysis) of a polypeptide at its peptide 
linkages. The amino acid residues described herein are presumed to be in 
the "L" isomeric form. Residues in the "D" isomeric form, which are so- 
designated, can be substituted for any L-amino acid residue, as long as 
the desired functional property is retained by the polypeptide. NH 2 refers 
to the free amino group present at the amino terminus of a polypeptide. 
COOH refers to the free carboxy group present at the carboxyl terminus 
of a polypeptide. In keeping with standard polypeptide nomenclature 
described in J. Biol. Chem., 243:3552-59 (1969) and adopted at 37 
C.F.R. § § 1.821 - 1.822, abbreviations for amino acid residues are 
shown in the following Table: 

Table 1 



Table of Correspondence 



SYMBOL 




1 -Letter 


3-Letter 


AMINO ACID 


Y 


Tyr 


tyrosine 


G 


Gly 


glycine 


F 


Phe 


phenylalanine 
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SYMBOL 1 




M 


Met 1 


methionine 


i. A 


Ala 1 


alanine 


S 


Ser 


serine 


1 1 


lie | 


isoleucine 


1 L 


Leu 


leucine 


I T 


Thr J 


threonine 


1 V 


Val | 


valine 


P 


Pro 


proline 


1 K 


Lys 


, lysine 


1 H 


His 


histidine 


Q 


Gin 


glutamine ' 


1 E 


Glu 


glutamic acid 


Z 


Glx 


Glu and/or Gin 


1 W 


Trp 


tryptophan 


1 R 


Arg 


arginine 


1 D 


Asp 


aspartic acid 


1 N 


Asn 


asparagine 


1 B 


Asx 


Asn and/or Asp 


1 C 


Cys 


cysteine 


1 X 


Xaa 


Unknown or other 



It should be noted that all amino acid residue sequences 
represented herein by formulae have a left to right orientation in the 
conventional direction of amino-terminus to carboxyl-terminus. In 
addition, the phrase "amino acid residue" is broadly defined to include the 
amino acids listed in the Table of Correspondence and modified and 
unusual amino acids, such as those referred to in 37 C.F.R. § § 1.821- 
1 .822, and incorporated herein by reference. Furthermore, it should be 
noted that a dash at the beginning or end of an amino acid residue 
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sequence indicates a peptide bond to a further sequence of one or more 
amino acid residues or to an amino-terminal group such as NH 2 or to a 
carboxyl-terminal group such as COOH. 

In a peptide or protein, suitable conservative substitutions of amino 
acids are known to those of skill in this art and may be made generally 
without altering the biological activity of the resulting molecule. Those of 
skill in this art recognize that, in general, single amino acid substitutions 
in non-essential regions of a polypeptide do not substantially alter 
biological activity (see, e.g., Watson et al. Molecular Biology of the Gene, 
4th Edition, 1987, The Benjamin/Cummings Pub. co., p. 224). 

Such substitutions are preferably made in accordance with those 
set forth in TABLE 2 as follows: 



TABLE 2 



Original residue 


Conservative substitution 


Ala (A) 


Gly; Ser 


Arg (R) 


Lys 


Asn (N) 


Gin; His 


Cys (C) 


Ser 


Gin (Q) 


Asn 


Glu (E) 


Asp 


Gly (G) 


Ala; Pro 


His (H) 


Asn; Gin 


He (I) 


Leu; Val 


Leu (L) 


lie; Val 


Lys (K) 


Arg; Gin; Glu 


Met (M) 


Leu; Tyr; lie 


Phe (F) 


Met; Leu; Tyr 


Pro (P) 


Ala; Gly 


Ser (S) 


Thr 


Thr (T) 


Ser 


Trp (W) 


Tyr 


Tyr (Y) 


Trp; Phe 


Val (V) 


lie; Leu 



Other substitutions are also permissible and may be determined 
empirically or in accord with known conservative substitutions. 

As used herein, nucleic acids include DNA, RNA and analogs 
thereof, including protein nucleic acids (PNA) and mixture thereof. 
Nucleic acids can be sjngle or double stranded. When referring to probes 
or primers, optionally labeled, with a detectable label, such as a 
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fluorescent or radiolabel, single-stranded molecules are contemplated. 
Such molecules are typically of a length such that they are statistically 
unique of low copy number (typically less than 5, preferably less than 3) 
for probing or priming a library. Generally a probe or primer contains at 
least 14, 16 or 30 contiguous of sequence complementary to or identical 
a gene of interest. Probes and primers can be 10, 14, 16, 20, 30, 50, 
100 or more nucleic acid bases long. 

As used herein, by homologous means about greater than 25% 
nucleic acid sequence identity, preferably 25% 40%, 60%, 80%, 90% or 
95%. The intended percentage will be specified. The terms "homology" 
and "identity" are often used interchangeably. In general, sequences are 
aligned so that the highest order match is obtained (see, e.g. : 
Computational Molecular Biology, Lesk, A.M., ed., Oxford University 
Press, New York, 1988; Biocomputing: Informatics and Genome Projects, 
Smith, D.W., ed., Academic Press, New York, 1993; Computer Analysis 
of Sequence Data, Part I, Griffin, A.M., and Griffin, H.G., eds., Humana 
Press, New Jersey, 1 994; Sequence Analysis in Molecular Biology, von 
Heinje, G., Academic Press, 1987; and Sequence Analysis Primer, 
Gribskov, M. and Devereux, J., eds., M Stockton Press, New York, 1991; 
Carillo et at. (1988) SI AM J Applied Math 45:1073). By sequence 
identity, the number of conserved amino acids are determined by standard 
alignment algorithms programs, and are used with default gap penalties 
established by each supplier. Substantially homologous nucleic acid 
molecules would hybridize typically at moderate stringency or at high 
stringency all along the length of the nucleic acid of interest. Also 
contemplated are nucleic acid molecules that contain degenerate codons 
in place of codons in the hybridizing nucleic acid molecule. 

As used herein, a nucleic acid homolog refers to a nucleic acid that 
includes a preselected conserved nucleotide sequence, such as a 
sequence encoding a therapeutic polypeptide. By the term "substantially 
homologous" is meant having at least 80%, preferably at least 90%, 



most preferably at least 95% homology therewith or a less percentage of 
homology or identity and conserved biological activity or function. 

The terms "homology" and "identity" are often used 
interchangeably. In this regard, percent homology or identity may be 
5 determined, for example, by comparing sequence information using a GAP 
computer program. The GAP program uses the alignment method of 
Needleman and Wunsch {J. Mol. Biol. 48:443 (1970), as revised by Smith 
and Waterman (Adv. Appl. Math. 2:482 (1981). Briefly, the GAP program 
defines similarity as the number of aligned symbols (i.e., nucleotides or 
0 amino acids) which are similar, divided by the total number of symbols in 
the shorter of the two sequences. The preferred default parameters for 
the GAP program may include: (1) a unitary comparison matrix 
(containing a value of 1 for identities and 0 for non-identities) and the 
weighted comparison matrix of Gribskov and Burgess, Nucl. Acids Res. 
5 14:6745 (1986), as described by Schwartz and Dayhoff, eds., ATLAS OF 
PROTEIN SEQUENCE AND STRUCTURE National Biomedical Research 
Foundation, pp. 353-358 (1979); (2) a penalty of 3.0 for each gap and an 
additional 0.10 penalty for each symbol in each gap; and (3) no penalty 
for end gaps. 

Whether any two nucleic acid molecules have nucleotide sequences 
that are, for example, at least 80%, 85%, 90%, 95%, 96%, 97%, 98% 
or 99% /'identical" can be determined using known computer algorithms 
such as the "FAST A" program, using for example, the default parameters 
as in Pearson and Lipman, Proc. Natl. Acad. Sci. USA 55:2444 (1988). 
Alternatively the BLAST function of the National Center for Biotechnology 
Information database may be used to determine identity 

In general, sequences are aligned so that the highest order match 
is obtained. "Identity" per se has an art-recognized meaning and can be 
calculated using published techniques. (See, e.g.: Computational 
Molecular Biology, Lesk, A.M., ed., Oxford University Press, New York, 
1988; Biocomputing: Informatics and Genome Projects, Smith, D.W., ed., 
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Academic Press, New York, 1993; Computer Analysis of Sequence Data, 
Part I, Griffin, A.M., and Griffin, H.G., eds., Humana Press, New Jersey, 
1994; Sequence Analysis in Molecular Biology, von Heinje, G., Academic 
Press, 1987; and Sequence Analysis Primer, Gribskov, M. and Devereux, 
5 J., eds., M Stockton Press, New York, 1991). While there exist a number 
of methods to measure identity between two polynucleotide or 
polypeptide sequences, the term "identity" is well known to skilled 
artisans (Carillo, H. & Lipton, D„ SI AM J Applied Math 45:1073 (1988)). 
Methods commonly employed to determine identity or similarity between 
0 two sequences include, but are not limited to, those disclosed in Guide to 
Huge Computers, Martin J. Bishop, ed., Academic Press, San Diego, 
1994, and Carillo, H. & Lipton, D., S/AM J Applied Math 45:1073 
(1988). Methods to determine identity and similarity are codified in 
computer programs. Preferred computer program methods to determine 
5 identity and similarity between two sequences include, but are not limited 
to, GCG program package (Devereux, J., et al., Nucleic Acids Research 
12(0:387 (1984)), BLASTP, BLASTN, FASTA (Atschul, S.F., et al., J 
Molec Biol 275:403 (1990)), and CLUSTALW. For sequences displaying 
a relatively high degree of homology, alignment can be effected manually 
by simply lining up the sequences and manually or visually matching the 
conserved portions. 

Therefore, as used herein, the term "identity" represents a 
comparison between a test and a reference polypeptide or polynucleotide. 
For example, a test polypeptide may be defined as any polypeptide that 
is 90% or more identical to a reference polypeptide. 

For the alignments presented herein (see, Fig. 5) for the AAV 
serotype, the CLUSTALW program was employed with parameters set as 
follows: scoring matrix BLOSUM, gap open 10, gap extend 0.1, gap 
distance 40% and transitions/transversions 0.5; specific residue penalties 
for hydrophobic amino acids (DEGKNPQRS), distance between gaps for 
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which the penalties are augmented was 8, and gaps of extremities 
penalized less than internal gaps. 

As used herein, a "corresponding" position on a protein, such as 
the AAV rep protein, refers to an amino acid position based upon 
alignment to maximize sequence identity. For AAV Rep proteins an 
alignment of the Rep 78 protein from AAV-2 and the corresponding 
protein from other AAV serotypes (AAV-1, AAV-6, AAV-3, AAV-3B, 
AAV-4, AAV-2 and AAV-5) is shown in Figure 5. The "hit" positions are 
shown in italics. 

As used herein, the term at least "90% identical to" refers to 
percent identities from 90 to 100% relative to the reference polypeptides. 
Identity at a level of 90% or more is indicative of the fact that, assuming 
for exemplification purposes a test and reference polynucleotide length of 
100 amino acids are compared. No more than 10% (i.e., 10 out of 100) 
amino acids in the test polypeptide differs from that of the reference 
polypeptides. Similar comparisons may be made between a test and 
reference polynucleotides. Such differences may be represented as point 
mutations randomly distributed over the entire length of an amino acid 
sequence or they may be clustered in one or more locations of varying 
length up to the maximum allowable, e.g. 10/100 amino acid difference 
(approximately 90% identity). Differences are defined as nucleic acid or 
amino acid substitutions, or deletions. 

As used herein, it is also understood that the terms substantially 
identical or similar varies with the context as understood by those skilled 
in the relevant art. 

As used herein, genetic therapy involves the transfer of 
heterologous nucleic acids to the certain cells, target cells, of a mammal, 
particularly a human, with a disorder or conditions for which such therapy 
is sought. The nucleic acid, such as DNA, is introduced into the selected 
target cells in a manner such that the heterologous nucleic acid, such as 
DNA, is expressed and a therapeutic product encoded thereby is 
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produced. Alternatively, the heterologous nucleic acid, such as DNA, 
may in some manner mediate expression of DNA that encodes the 
therapeutic product, or it may encode a product, such as a peptide or 
RNA that in some manner mediates, directly or indirectly, expression of a 
therapeutic product. Genetic therapy may also be used to deliver nucleic 
acid encoding a gene product that replaces a defective gene or 
supplements a gene product produced by the mammal or the cell in which 
it is introduced. The introduced nucleic acid may encode a therapeutic 
compound, such as a growth factor inhibitor thereof, or a tumor necrosis 
factor or inhibitor thereof, such as a receptor therefor, that is not normally 
produced in the mammalian host or that is not produced in therapeutically 
effective amounts or at a therapeutically useful time. The heterologous 
nucleic acid, such as DNA, encoding the therapeutic product may be 
modified prior to introduction into the cells of the afflicted host in order to 
enhance or otherwise alter the product or expression thereof. Genetic 
therapy may also involve delivery of an inhibitor or repressor or other 
modulator of gene expression. 

As used herein, heterologous or foreign nucleic acid, such as DNA 
and RNA, are used interchangeably and refer to DNA or RNA that does 
not occur naturally as part of the genome in which it is present or which 
is found in a location or locations in the genome that differ from that in 
which it occurs in nature. Heterologous nucleic acid is generally not 
endogenous to the cell into which it is introduced, but has been obtained 
from another cell or prepared synthetically. Generally, although not 
necessarily, such nucleic acid encodes RNA and proteins that are not 
normally produced by the cell in which it is expressed. Any DNA or RNA 
that one of skill in the art would recognize or consider as heterologous or 
foreign to the cell in which it is expressed is herein encompassed by 
heterologous DNA. Heterologous DNA and RNA may also encode RNA or 
proteins that mediate or alter expression of endogenous DNA by affecting 
transcription, translation, or other regulatable biochemical processes. 
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Examples of heterologous nucleic acid include, but are not limited to, 
nucleic acid that encodes traceable marker proteins, such as a protein 
that confers drug resistance, nucleic acid that encodes therapeutically 
effective substances, such as anti-cancer agents, enzymes and hormones, 
and DNA that encodes other types of proteins, such as antibodies. 

Hence, herein heterologous DNA or foreign DNA, includes a DNA 
molecule not present in the exact orientation and position as the 
counterpart DNA molecule found in the genome. It may also refer to a 
DNA molecule from another organism or species (i.e., exogenous). 

As used herein, a therapeutically effective product introduced by 
genetic therapy is a product that is encoded by heterologous nucleic acid, 
typically DNA, that, upon introduction of the nucleic acid into a host, a 
product is expressed that ameliorates or eliminates the symptoms, 
manifestations of an inherited or acquired disease or that cures the 
disease. 

As used herein, isolated with reference to a nucleic acid molecule 
or polypeptide or other biomolecule means that the nucleic acid or 
polypeptide has separated from the genetic environment from which the 
polypeptide or nucleic acid were obtained. It may also mean altered from 
the natural state. For example, a polynucleotide or a polypeptide naturally 
present in a living animal is not "isolated," but the same polynucleotide or 
polypeptide separated from the coexisting materials of its natural state is 
"isolated", as the term is employed herein. Thus, a polypeptide or 
polynucleotide produced and/or contained within a recombinant host cell 
is considered isolated. Also intended as an "isolated polypeptide" or an 
"isolated polynucleotide" are polypeptides or polynucleotides that have 
been purified, partially or substantially, from a recombinant host cell or 
from a native source. For example, a recombinantly produced version of 
a compounds can be substantially purified by the one-step method 
described in Smith and Johnson, Gene 57:31-40 (1988). The terms 
isolated and purified are sometimes used interchangeably. 
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Thus, by "isolated" is meant that the nucleic is free of the coding 
sequences of those genes that, in the naturally-occurring genome of the 
organism (if any) immediately flank the gene encoding the nucleic acid of 
interest. Isolated DNA may be single-stranded or double-stranded, and 
may be genomic DNA, cDNA, recombinant hybrid DNA, or synthetic 
DNA. It may be identical to a native DNA sequence, or may differ from 
such sequence by the deletion, addition, or substitution of one or more 
nucleotides. 

Isolated or purified as it refers to preparations made from biological 
cells or hosts means any cell extract containing the indicated DNA or 
protein including a crude extract of the DNA or protein of interest. For 
example, in the case of a protein, a purified preparation can be obtained 
following an individual technique or a series of preparative or biochemical 
techniques and the DNA or protein of interest can be present at various 
degrees of purity in these preparations. The procedures may include for 
example, but are not limited to, ammonium sulfate fractionation, gel 
filtration, ion exchange change chromatography, affinity chromatography, 
density gradient centrifugation and electrophoresis. 

A preparation of DNA or protein that is "substantially pure" or 
"isolated" should be understood to mean a preparation free from naturally 
occurring materials with which such DNA or protein is normally 
associated in nature. "Essentially pure" should be understood to mean a 
"highly" purified preparation that contains at least 95% of the DNA or 
protein of interest. 

A cell extract that contains the DNA or protein of interest should be 
understood to mean a homogenate preparation or cell-free preparation 
obtained from cells that express the protein or contain the DNA of 
interest. The term "cell extract" is intended to include culture media, 
especially spent culture media from which the cells have been removed. 

As used herein, receptor refers to a biologically active molecule 
that specifically binds to (or with) other molecules. The term "receptor 



-23- 



protein" may be used to more specifically indicate the proteinaceous 
nature of a specific receptor. 

As used herein, recombinant refers to any progeny formed as the 
result of genetic engineering. 

As used herein, a promoter region refers to the portion of DNA of a 
gene that controls transcription of the DNA to which it is operatively 
linked. The promoter region includes specific sequences of DNA that are 
sufficient for RNA polymerase recognition, binding and transcription 
initiation. This portion of the promoter region is referred to as the 
promoter. In addition, the promoter region includes sequences that 
modulate this recognition, binding and transcription initiation activity of 
the RNA polymerase. These sequences may be cis acting or may be 
responsive to trans acting factors. Promoters, depending upon the nature 
of the regulation, may be constitutive or regulated. 

As used herein, the phrase "operatively linked" generally means the 
sequences or segments have been covalently joined into one piece of 
DNA, whether in single or double stranded form, whereby control or 
regulatory sequences on one segment control or permit expression or 
replication or other such control of other segments. The two segments 
are not necessarily contiguous. For gene expression a DNA sequence and 
a regulatory sequence(s) are connected in such a way to control or permit 
gene expression when the appropriate molecular, e.g., transcriptional 
activator proteins, are bound to the regulatory sequence(s). 

As used herein, production by recombinant means by using 
recombinant DNA methods means the use of the well known methods of 
molecular biology for expressing proteins encoded by cloned DNA, 
including cloning expression of genes and methods, such as gene 
shuffling and phage display with screening for desired specificities. 

As used herein, a splice variant refers to a variant produced by 
differential processing of a primary transcript of genomic DNA that results 
in more than one type of mRNA. 
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As used herein, a composition refers to any mixture of two or more 
products or compounds. It may be a solution, a suspension, liquid, 
powder, a paste, aqueous, non-aqueous or any combination thereof. 

As used herein, a combination refers to any association between 
two or more items. 

As used herein, substantially identical to a product means 
sufficiently similar so that the property of interest is sufficiently 
unchanged so that the substantially identical product can be used in place 
of the product. 

As used herein, the term "vector" refers to a nucleic acid molecule 
capable of transporting another nucleic acid to which it has been linked. 
One type of preferred vector is an episome, i.e., a nucleic acid capable of 
extra-chromosomal replication. Preferred vectors are those capable of 
autonomous replication and/or expression of nucleic acids to which they 
are linked. Vectors capable of directing the expression of genes to which 
they are operatively linked are referred to herein as "expression vectors". 
In general, expression vectors of utility in recombinant DNA techniques 
are often in the form of "plasmids" which refer generally to circular 
double stranded DNA loops which, in their vector form are not bound to 
the chromosome. "Plasmid" and "vector" are used interchangeably as the 
plasmid is the most commonly used form of vector. Other such other 
forms of expression vectors that serve equivalent functions and that 
become known in the art subsequently hereto. 

As used herein, vector is also used interchangeable with "virus 
vector" or "viral vector". In this case, which will be clear from the 
context, the "vector" is not self-replicating. Viral vectors are engineered 
viruses that are operatively linked to exogenous genes to transfer (as 
vehicles or shuttles) the exogenous genes into cells. 

As used herein, transduction refers to the process of gene transfer 
and expression into mammalian and other cells mediated by viruses. 
Transfection refers to the process when mediated by plasmids. 
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As used herein, "polymorphism" refers to the coexistence of more 
than one form of a gene or portion thereof. A portion of a gene of which 
there are at least two different forms, i.e., two different nucleotide 
sequences, is referred to as a "polymorphic region of a gene". A 
polymorphic region can be a single nucleotide, referred to as a single 
nucleotide polymorphism (SNP), the identity of which differs in different 
alleles. A polymorphic region can also be several nucleotides in length. 

As used herein, "polymorphic gene" refers to a gene having at least 
one polymorphic region. 

As used herein, "allele", which is used interchangeably herein with 
"allelic variant" refers to alternative forms of a gene or portions thereof. 
Alleles occupy the same locus or position on homologous chromosomes. 
When a subject has two identical alleles of a gene, the subject is said to 
be homozygous for the gene or allele. When a subject has two different 
alleles of a gene, the subject is said to be heterozygous for the gene. 
Alleles of a specific gene can differ from each other in a single nucleotide, 
or several nucleotides, and can include substitutions, deletions, and 
insertions of nucleotides. An allele of a gene can also be a form of a gene 
containing a mutation. 

As used herein, the term "gene" or "recombinant gene" refers to a 
nucleic acid molecule comprising an open reading frame and including at 
least one exon and (optionally) an intron sequence. A gene can be either 
RNA or DNA. Genes may include regions preceding and following the 
coding region (leader and trailer). 

As used herein, "intron" refers to a DNA sequence present in a 
given gene which is spliced out during mRNA maturation. 

As used herein, "nucleotide sequence complementary to the 
nucleotide sequence set forth in SEQ ID NO: x" refers to the nucleotide 
sequence of the complementary strand of a nucleic acid strand having 
SEQ ID NO: x. The term "complementary strand" is used herein 
interchangeably with the term "complement". The complement of a 
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nucleic acid strand can be the complement of a coding strand or the 
complement of a non-coding strand. When referring to double stranded 
nucleic acids, the complement of a nucleic acid having SEQ ID NO: x 
refers to the complementary strand of the strand having SEQ ID NO: x or 
5 to any nucleic acid having the nucleotide sequence of the complementary 
strand of SEQ ID NO: x. When referring to a single stranded nucleic acid 
having the nucleotide sequence SEQ ID NO: x, the complement of this 
nucleic acid is a nucleic acid having a nucleotide sequence which is 
complementary to that of SEQ ID NO: x. 
0 As used herein, the term "coding sequence" refers to that portion 

of a gene that encodes an amino acid sequence of a protein. 

As used herein, the term "sense strand" refers to that strand of a 
double-stranded nucleic acid molecule that has the sequence of the 
mRNA that encodes the amino acid sequence encoded by the double- 
5 stranded nucleic acid molecule. 

As used herein, the term "antisense strand" refers to that strand of 
a double-stranded nucleic acid molecule that is the complement of the 
sequence of the mRNA that encodes the amino acid sequence encoded 
by the double-stranded nucleic acid molecule. 

As used herein, an array refers to a collection of elements, such as 
nucleic acid molecules, containing three or more members. An 
addressable array is one in which the members of the array are 
identifiable, typically by position on a solid phase support or by virtue of 
an identifiable or detectable label, such as by color, fluorescence, 
electronic signal (i.e. radiofrequency (RF), microwave or other frequency 
that does not substantially alter the interaction of the molecules of 
interest), bar code or other symbology, chemical or other such label. 
Hence, in general the members of the array are immobilized to discrete 
identifiable loci on the surface of a solid phase or directly or indirectly 
linked to or otherwise associated with the identifiable label, such as 
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affixed to a microsphere or other particulate support (herein referred to as 
beads) and suspended in solution or spread out on a surface. 

As used herein, a support (also referred to as a matrix support, a 
matrix, an insoluble support or solid support) refers to any solid or 
semisolid or insoluble support to which a molecule of interest, typically a 
biological molecule, organic molecule or biospecific ligand is linked or 
contacted. Such materials include any materials that are used as affinity 
matrices or supports for chemical and biological molecule syntheses and 
analyses, such as, but are not limited to: polystyrene, polycarbonate, 
polypropylene, nylon, glass, dextran, chitin, sand, pumice, agarose, 
polysaccharides, dendrimers, buckyballs, polyacrylamide, silicon, rubber, 
and other materials used as supports for solid phase syntheses, affinity 
separations and purifications, hybridization reactions, immunoassays and 
other such applications. The matrix herein can be particulate or can be 
in the form of a continuous surface, such as a microtiter dish or well, a 
glass slide, a silicon chip, a nitrocellulose sheet, nylon mesh, or other 
such materials. When particulate, typically the particles have at least one 
dimension in the 5-10 mm range or smaller. Such particles, referred 
collectively herein as "beads", are often, but not necessarily, spherical. 
Such reference, however, does not constrain the geometry of the matrix, 
which may be any shape, including random shapes, needles, fibers, and 
elongated. Roughly spherical "beads", particularly microspheres that can 
be used in the liquid phase, are also contemplated. The "beads" may 
include additional components, such as magnetic or paramagnetic 
particles (see, e.g.,, Dyna beads (Dynal, Oslo, Norway)) for separation 
using magnets, as long as the additional components do not interfere with 
the methods and analyses herein. 

As used herein, matrix or support particles refers to matrix 
materials that are in the form of discrete particles. The particles have any 
shape and dimensions, but typically have at least one dimension that is 
100 mm or less, 50 mm or less, 10 mm or less, 1 mm or less, 100 jjm or 
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less, 50 pm or less and typically have a size that is 100 mm 3 or less, 50 
mm 3 or less, 10 mm 3 or less, and 1 mm 3 or less, 100 //m 3 or less and may 
be order of cubic microns. Such particles are collectively called "beads." 

As used herein, the abbreviations for any protective groups, amino 
acids and other compounds, are, unless indicated otherwise, in accord 
with their common usage, recognized abbreviations, or the IUPAC-IUB 
Commission on Biochemical Nomenclature (see, (1972) Biochem. 
7 7:942-944). 

B. High Throughput Process 

Provided herein are high throughput process for the generation of 
and identification of proteins that exhibit desired phenotypes. The 
processes, include methods that are particularly adapted for high 
throughput protocols, which require accurate methods for identifying 
modified proteins. 

A general directed evolution process includes the following steps: 

1 . Generation of diversity at the nucleic acid level, on the gene 
to be 'evolved' 

2. phenotypic characterization of the gene variants generated; 

and 

3. identification of optimized gene variants. 

The processes provided herein effect these steps that can be 
performed in a high throughput format (see, FIGURES 1) that is optionally 
automated. A distinguishing characteristic of the processes provided 
herein, is that each candidate nucleic acid molecule is separately 
generated and screened. In an automated process at least some of the 
steps are performed without human intervention and are generally 
controlled by software. Most, if not all steps, are performed in 
addressable formats, such as at discrete locations in or on solid supports, 
such as microtiter plates or in other addressable formats, such as linked 
to coded supports. The supports can be electronically, physically, 
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chemically or otherwise identifiable, such as by an identifiable symbology, 
including n bar code, or can be color coded. 

1 . Gen ration of Diversity using a semi-rational approach 

A semi-rational approach to creating diversity or evolving genes is 
provided herein. The goal is to create diversity but to decrease the 
number of molecules to be screened. By reducing the numbers, the 
molecules can be screened in high throughput format molecule-by- 
molecule (or groups thereof). 

Generation of diversity at the nucleic level, in principle, can be 
accomplished by a number of diverse technologies like mutagenesis 
(either site-directed or random), recombination, shuffling and de-novo 
synthesis. These different technologies differ in the degree of diversity 
they generate as well as in the minimal length of the unitary change they 
can introduce (from single base to large domains). The outcome of step 1 
is a collection of diverse, although highly related, molecules that 
constitutes what is called a 'library'. 

This step is crucial, since it provides the initial conditions for the 
entire process and is determinative of the outcome. The chances of 
finding an optimized gene version in a library is a function of the total 
diversity present in the library. In addition, the type of diversity 
introduced (such as, but not limited to, single point mutations, multiple 
point mutations, scarce small rearrangements, recombination of large 
domains, multiple shuffling) condition the outcome, particularly with 
respect to the generation of new variants compared to the original gene, 
and the probability that the new variants, not only exhibit the "evolved" 
function or property, but also work in their natural biological networks 
where they are expected to act by interacting, recognizing, and/or being 
recognized, by a large panoply of other proteins and other molecules. 
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Rapid disc very f protein variants at th amino acid level by 
rational mutagen sis (aa-scan) 

A method, referred to herein as an amino-acid scan method for 
directed evolution, is provided herein for generating protein variants. This 
5 method can be performed on an entire protein or selected domains 
thereof, or can be used to identify benchmark sequences, such as 
functional domains, and, for example, recombine them as exchangeable 
units or restrict the diversity to limited or specific regions of the protein. 
M: Not only can this method be used with the processes provided herein, but 

10 it also has applications for any methods that use such variants or require 
generation of such variants, such as, but not limited to, searches for 
consensus sequences and homology regions that are used in functional 
genomics, functional proteomics; comparative modeling in protein 
crystallography and protein modeling; searches for natural diversity, 
15 (e.g., directed evolution methods in 6,171,820, 6,238,884, 6,174,673, 
6,057,103, 6,001,574, 5,763,239, ); exon- or family-shuffling-based 
diversity (e.g., directed evolution using gene shuffling (see, e.g., U.S. 
Patent Nos. 6,096,548, 6,117,679, 6,165,793, 6,180,406, 6,132,970); 
the optimization of only the CDRs regions (e.g., directed evolution of 
20 antibodies see., e.g., U.S. Patent Nos. 5,723,323, 6,258,530, 5770,434, 
5,862,514) and other methods (see, e.g., U.S. Patent Nos. 5,837,500, 
5,571,698, 6,156,509). 

The amino-acid scanning-based method provided herein has 
advantages that prior methods do not have. For example, prior methods 
25 are based upon the underlying assumption that there are parts of the 
molecule (gene or protein) that are sufficiently adapted to perform their 
respective function, and further changes are not advantageous. Such 
methods do not look at total potential plasticity of a given molecule, but 
at the plasticity still permitted while keeping some basic functions in 
30 place. By choosing this route, however, additional potential variation is 
missed. The potential in the intrinsic plasticity of those regions that are 
presumed 'preserved' is lost. For instance, methods (e.g., those in U.S. 
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Patent Nos. 6,171,820, 6,238,884, 6,174,673, 6,057,103, 6,001,574, 
5,763,239) that use natural diversity can miss the potential plasticity 
within those regions that are naturally 'conserved', i.e. there where there 
is no natural diversity. Methods that rely on exon- or family-shuffling- 
based diversity can miss the potential plasticity within regions contained 
in the shuffled fragments, i.e. within the fragments exchanged as a block. 

The method provided herein in contrast is sufficiently flexible to 
create mutants at a variety of levels, including at the single amino acid 
level; i.e. the method can generate mutants that differ from each other at 
a single amino and not at a larger block level. The challenge solved by 
the method herein is to generate diversity at the single amino acid level, 
without moving too close to a pure 'random' approach, which results in 
an intractable number of mutants. 

The method provided herein is based on the premise that there are 
single amino acids or small blocks of sequence of amino acids that are 
either (1) directly involved in the activities that the methods 'evolve' 
(these amino acids would be at or close to the 'active sites' of the 
protein); or (2) directly involved in maintaining within the protein the intra- 
molecular environment that allows the active site(s) to stay active. 

Potential plasticity at the amino acid level can be exploited if amino 
acids or blocks of amino acids directly involved in the active sites for the 
activity to be evolved are known. Often they are not known. The 
problem that is solved herein, however, is how to exploit the potential 
plasticity at the amino acid level when nothing is known about the 
structure of the protein in question or about the position of its single or 
numerous active sites. 

The technology referred to herein as amino acid-scanning has been 
used to precisely identify those amino acids directly involved in the active 
sites of some enzymes and receptors (see, e.g., Becl-Sickinger et al. 
(1994) Eur. J. Biochem. 223:947-958; Gibbs et al. (1991)J. Biol. Chem. 
256:8923-8931; Matsushita et al. (2000) J. Biol. Chem. 275:11044- 
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1 1049) but has not been employed for directed evolution or for the 

generation of diversity. The amino acid scan as practiced in the prior art 

is used to produce a set of protein mutants, often within the region 

* 

suspected to contain the active site(s), such that in each individual 
5 mutant a selected residue, such as Ala, replaces a different amino acid. 
Ala or other neutral amino acids generally, although not necessarily, is 
selected as a replacement amino acid since, except in instances in which 
the replaced amino acid is directly involved in an active site, it should 
have a neutral effect on the protein activity and not disturb the native 

10 secondary structure of the protein. In instances in which the replaced 

amino acid is directly involved in an active site the activity of that site will 
be lost or altered. Amino-acid scanning, such as Ala-scanning, has been 
successfully applied for the identification of active sites in a number of 
proteins, and has been performed in computer-based rational drug design 

15 methods. Other amino acids, particularly amino acids that have a neutral 
effect, such as Glycine, can also be used. 

The amino acid scanning method is employed herein for the 
generation of the mutant proteins for screening for identification of sites 
or loci in a target protein or regions in a protein that alter a selected 

20 activity. In performing this method, the amino acids are each replaced, 
one-by-one along the full-length of the protein, or one-by-one in pre- 
selected domains, such as domains that possess a desired activity or 
exhibit a particular function. Once sites of interest are identified other 
methods for generating diversity from the resulting molecules can be 

25 employed or the further steps of the method provided herein can be 
performed. 

The method includes the following steps: 

(1) Identification of the active site(s) on the full 
length protein sequence. In one embodiment a full-length amino acid- 

30 scan, typically, although not necessarily, an Ala-scan, or the identification 
and positioning of the active site(s) on proteins of either known or 



-33- 



unknown function. For purposes herein, an active site is not necessarily 
the natural active sites involved in the natural activity of a target protein, 
but those amino acids involved in the activities of the proteins under 
'directed evolution' with the purpose of either gain, improvement or loss 
of function. 

The whole process of the 'identification of the active site(s) on the 
full length protein sequence requires the following sub-steps: 

a. Generation of a mutant library (on the 
gene to be evolved) in which each individual mutant contains a single 
mutation located at a different amino acid position and that includes a 
systematic replacement of the native amino acid by Ala or any other 
amino acid (always the same throughout the entire protein sequence); 

b. phenotypic characterization of the 
individual mutants, one-by-one and assessment of mutant protein activity; 

c. identification of those mutants that 
display an alteration, typically a decrease, in the selected protein activity, 
thus, indicating that amino acids directly involved in the active site(s) 
have been hit. The aa positions whose aa-scan mutations display an 
alteration, typically a loss or decrease, in activity are named HITS. 

The identification of the active site(s) (HITS) is thus, by this 
method, made in a completely unbiased manner. There are no 
assumptions about the specific structure of the protein in question nor 
any knowledge or assumptions about the active site(s). The results of the 
amino acid scan identify such sites. 

Once the active site(s) (the HITS) has(ve) been identified, those 
amino acids either at or surrounding the active sites, such as within 1, 2, 
3, ... 10, 20 or any selected regions, as the unitary elements of 
exchange and generate diversity either at or around one of those sites or 
as a combined diversity at several sites at a time can be assessed. 
This process includes the following steps: 
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a. Generation of a new mutant library (on the gene to be 
evolved) in which each individual mutant contains either single or multiple 
mutations located at (or surrounding) a specific active site (a HIT) position 
detected by the precedent aa-scan process. In the example these 
mutations include replacement, in each individual mutant, of the native 
amino acid located either at (or surrounding) the HIT position by one of all 
other possible amino acids, such that, in the library, and at (or 
surrounding) each HIT position the native amino acid has been replaced 
by all possible amino acids. 

b. Identification of those mutants that display an 
increase in protein activity, thus indicating that a new sequence at or 
surrounding an active site has been identified with higher activity 
compared to the native sequence. These optimized sequences are named 
LEADS. 

The process can be repeated as many times as desired, in search 
for new combinations of optimal amino acids at (or surrounding) the 
different HIT sites. Each time, the process includes the steps of 
generating of a new mutant library (of the gene to be evolved) in which 
each individual mutant contains either single or multiple mutations located 
at (or surrounding) a specific active site (a HIT) position; phenotypic 
characterization of the individual mutants, one-by-one and assessment of 
mutant protein activity; and identification of those mutants that display an 
increase in protein activity, thus indicating that a new sequence at or 
surrounding an active site has been identified with higher activity 
compared to the native sequence. These optimized sequences are again 
named LEADS (second generation LEADS). 

2. Phenotypic characterization of the gene variants 
This step requires the expression of the gene variants in order to 
allow them to manifest their respective phenotypes. Gene expression can 
be accomplished by different means: in vitro, in reconstituted systems or 
in vivo in cellular systems, including bacterial and eukaryotic cells. For all 
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exemplification purposes, reference is made to in vivo systems. Those 
of skill in the art can readily adapt these methods for in vitro systems, 
including those using biochemical assays. 

This step is a crucial step for several reasons: 

(a) Expression system and protein processing. 
Depending on the system used (either bacteria or eukaryotic cells), 

as well as on the specific gene to be 'evolved', the variant proteins may 
or may not be appropriately processed, especially when post-translational 
modifications are involved, and therefore be able or not to elicit their 
potential activity. Consequently, the expression system (bacteria vs. 
eukaryotic cells) has to be carefully chosen. 

(b) Standardization of the expression system. 

The technologies available for gene transfer and expression into 
either bacteria or eukaryotic (let's consider mammalian) cells widely vary 
in their intrinsic efficiencies. While it is very easy to efficiently transfer 
and express genes in bacteria by chemical/physical methods 
(transformation), that is not the case for mammalian cells, where the 
transformation (here called transfection) process is inefficient and 
unreliable, specially when reproducibility and robustness are necessary in 
miniature, large number- and small scale high throughput settings like 
those necessary to analyze gene variant libraries. Therefore, when 
transfection is used on mammalian cells, the specific activity measured 
for the individual variants in the library most probably does not accurately 
reflect the real specific activity of the molecules involved. As provided 
herein, transduction, the process of gene transfer and expression into 
mammalian and other cells mediated by viruses, overcomes the limitations 
of transfection. 
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(3) Characterization. 

A distinction must be made between the 'expression' of the gene 
variants and their 'phenotypic characterization'. The expression system 
(either bacteria or mammalian cells) is only the vehicle to convert the 
gene variants into protein variants. The phenotypic characterization is 
performed on the protein variants, and may have nothing to do, 
depending on the specific system under study, with the cellular system 
used to express the variants. The phenotypic characterization requires the 
use of specific assays (either biochemical (cell-free) or cell-based assays) 
in which the activity of the different cell mutants is challenged and 
assessed. In addition to the implications discussed below, these assays 
must be designed in such a way that they reflect the final environment in 
which the 'evolved' protein is expected to act. As an example, when 
optimizing an enzyme to be used in an artificial industrial setting, the 
assay should reproduce those conditions (temperature, pH, media 
composition...) of the real-life industrial reaction mixture, which may be 
relatively easy to do. When the final destination of the 'evolved' protein 
is a complex biological setting, such as the intracellular environment, the 
extra-cellular milieu (example: circulating proteins) or the structure of a 
virus, the necessary assay(s) may be quite difficult to setup. With a few 
exceptions, most of the work done so far on directed evolution has been 
made on simple enzymes for which all the necessary settings are 
relatively easy to implement. 

Methods for accurately titering viruses 

Much progress in gene therapy, genomics, biotechnology and, in 
general, biomedical sciences, depends on the ability to generate and 
analyze large numbers and small amounts of specific viruses. High 
throughput technologies are employed in disciplines such as functional 
genomics and gene therapy in which the use of viruses plays a key role 
for the efficient transfer and analysis of large collections of genes or 
libraries. Also, virus samples and biomedical samples containing viruses 
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are routinely analyzed in thousands of hospitals, health centers, 
academic labs and biotech setting. 

Furthermore in processes herein, accurate titration can be important 
in at least two steps in the process. After preparation of the viruses with 
the mutated variant, and prior to screening, it is necessary to know the 
concentration of titer of the viruses in the sample, so that results among 
the samples can be compared. The methods in this section designated 
Real Time Virus Titering (RTVT™) and (TREE) are advantageously used. 

The methods in this section are also used in data analysis when 
measuring the output signal. As described below, output signal can be 
assessed by a Hill analysis or a second order polynomial or other 
algorithm that describes the interaction of biological molecules in complex 
system. In addition, where the output signal is actually the number of 
viral particles or ip produced, the methods in this section RTVT and TREE 
are advantageously used. 

Prior art virus titration methods (RCA, dRA...),for determining the 
amount of virus present in a biological sample, are based on the 
assessment of some kind of output signal, such as cytopathic effect, lysis 
or plaques and cell fusion focuses, induced in a reporter cell following a 
fixed time after infection with varying concentrations of the sample 
containing the virus. The lowest concentration of the sample at which no 
signal can be measured is taken as the titer of the virus in that sample. 
These approaches are known as "serial dilution" or "limiting dilution" 
methods. In limiting dilution methods, one single virus concentration, 
measured at a given time end-point gives rise to a single measurement of 
the output signal. These methods tend to be destructive in that 
assessment destroys the reaction so that only a single measurement can 
be taken on a sample. 
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Real Time Virus Tit ring (RTVT™) 

When a virus infects a cell, the infected cell undergoes a number 
of changes that can be followed over time and quantified. Such changes 
are designated herein as the "output signal". The cell reports an output 
signal in response to the infection and, therefore, it is named here a 
reporter cell. One such output signal, is, for example, the expression of 
the genes carried by the virus (whether they are viral genes or 
exogenous genes (transgenes)). The output signal (for instance the level 
of expression of those genes) develops over time and depends, mainly, 
on two factors: i) the time point ("t") at which its level is measured after 
infection and ii) the amount of virus infecting the cell; i.e. the 
concentration of the virus preparation used to infect the cell ("s"). 
The output signal, at a given time point after infection, will be higher for 
higher concentrations of the virus infecting the reporter cells; and for any 
given concentration of virus, the output signal increases with time after 
infection until it generally reaches a plateau level. 

Real Time Virus Titering (RTVT) published as International PCT 
application No. WO 01/186291, which is based on PCT/FR01/01366 and 
EXAMPLES below) uses non-destructive methods for the assessment of 
output signal. Real Time Virus Titering is a viral titration method based on 
the kinetic analysis of the development of the output signal in • 
virus-infected cells, tested at a single concentration of virus or biological 
sample. Instead of fixing the time point after infection and varying the 
concentration of the sample as is done in limiting dilution methods, in the 
Real Time Virus Titerng RTVT™ method, a fixed concentration of virus is 
used and the generation of a signal over time is assessed. Hence the 
signal is measured as a function of time at a single virus concentration. In 
this situation, a single virus sample (concentration), whose output signal 
is measured at a number of time points, can give rise to as many 
measurements of the output signal as needed, and, eventually to a 
continuous, over time, assessment of the signal in real time. 
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Real Time Virus Titerng RTVT can be advantageously used in high 
throughput methods in which large numbers of biological samples are 
analyzed at the same time. This is the case, for instance, when titering 
viruses in a virus library. Limiting dilution methods rely on the output 
signal generated by a number of dilutions of each individual sample. If, for 
example, 10 dilutions (or experimental points) of each virus are used for a 
titration using a limiting dilution method, the analysis of a library 
containing 10,000 viruses require analysis of 10 5 (i.e., 10 x 10,000) 
experimental points. The Real Time Virus Titerng RTVT™ method requires 
only one dilution per sample, thereby requiring 10-fold fewer experimental 
points than a limiting dilution method. For a Real Time Virus Titerng 
RTVT'" titering system, the time necessary for the output signal to 
reach a reference value (J3) is a direct function of the concentration of 
virus. Thus, t/3 can be used to directly determine the concentration of the 
virus. 

A limitation of the Real Time Virus Titerng RTVT 1 " limiting dilution 
titering method, however, is that not all the viruses (nor the genes carried 
by the viruses) generate a readily measured output signal that can be 
followed over time using non-destructive methods. 

Tagged Replication and expression enhancement 
(TREE) 

A method for titering designated Tagged Replication and Expression 
Enhancement Technology (TREE'") is provided herein. This system 
includes: i) a cell, ii) a reporter virus carrying a reporter gene, whose 
activity can be followed over time by a non- destructive method (i.e., 
fluorescence), iii) the virus (or virus library to be titered), herein referred 
to as the "titering virus". The elements are employed such that the virus 
to be titered interferes with any output signal generated by the reporter 
virus, leading to either decrease or increase in the amount of that signal. 
The higher the amount of virus to be titered, the higher is the interference 
with the reporter virus and output signal. In the absence of virus to be 
titered, the kinetics of the output signal generated by the reporter virus 
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are followed using the Real Time Virus Titemg RTVT titering method. In 
the presence of increasing amounts of the virus to be titered the output 
signal takes longer (or shorter) to develop as a function of the amount of 
virus to be titered. 

Using the TREE™ titering method, tfi, the time necessary for the 
output signal to reach a reference value (fi) is a direct function of the 
concentration of virus and, therefore, tfi can be used to determine the 
concentration of the virus to be titered. It is demonstrated herein (see the 
EXAMPLES) that when using the TREE system for titering, once an 
appropriate reference value {fi) is determined for the output signal 
generated from the reporter virus, the time tfi is a function of the 
concentration of the virus being titered (see Example). Therefore, the 
concentration (titer) of the virus to be titered, can be assessed by 
assessing the change induced in tfi by an aliquot of the virus to be titered. 
In a calibrated TREE titering assay, only one aliquot virus to be titered is 
needed to determine its titer, which is determined by measuring the shift 
in the tfi of the system. The only condition is that the virus to be titered 
must "interfere" (i.e., increase or decrease) the output signal of the 
reporter virus. 

A calibration curve representing tfi vs. the amount of virus to be 
titered is obtained using aliquots of a reference batch of virus of known 
titer (previously determined using any titering procedure). The calibration 
curve can then be used to determine the amount of virus in a sample of 
unknown titer, based on the change caused by an aliquot of the sample 
on the tfi of the system and the corresponding titer read from the 
calibration curve. 

3. Identification of gene variants 

There are at least two considerations in this step: 
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(a) S I ction vs. sere ning. 

Depending on the specific protein involved, and under certain and 
very specific assay conditions, those variants that have been 'evolved' 
may elicit a selective advantage over the native version. This situation 
represents the most simple case: the cells (bacteria or mammalian) 
expressing the library of protein variants, as a pool or mixture, can be 
simply exposed to the selective conditions which by themselves will allow 
to put in evidence the best optimized variants. This situation is however 
very rare and difficult to achieve. It can be hardly believed that for any 
protein that one may want to optimize, a suitable 'selective' assay could 
be set up. For the vast majority of the cases, selection will not be 
possible. Therefore pools of molecules cannot be used, because the 
specific readouts of the assays could not be attributed to individual 
variants. When the simplistic selection approach is not possible, then two 
things are absolutely needed: (a) a 'one-by-one' approach, i.e. each 
individual variant must be physically separated from the others and its 
activity tested independently; (b) an accurate and quantitative analysis 
that can distinguish slight differences in activity among the different 
variants along a wide range of performance values. 

(b) Accurate quantitative analytics 

When selection is not possible, the optimized variants must be 
distinguished from the native variant otherwise. The different degrees of 
optimization among the different variants in the library should, in addition, 
be distinguished if those variants showing the highest optimization level 
are to be identified. A powerful quantitative analytical protocol is then 
mandatory. These analytics should be able to attribute quantitative 
features (on the activity tested in the specific assay) to each of the 
variants tested and to rank them according to their individual 
performance. This requires, in addition, that each variant in the library is 
assayed individually; the use of pools or mixtures of molecules would 
hamper the ability to identify the right variants. 
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For such analysis, the output signal can be assessed by a Hill 
analysis (see Examples and (published International PCT application No. 
WO 01/44809 based on PCT n° PCT/FROO/03503),.) or a second order 
polynomial (see, Examples and (Drittanti eta/. (2000) Gene Then 7: 924- 
5 929)) or other algorithm that describes the interaction of biological 

molecules in complex system, such as the interaction between cells and 
biological agents. In addition, where the output signal is actually the 
number of viral particles or ip produced, the methods in this designated 
Real Time Virus Titering (RTVT™) and Tagged Replication and expression 
0 enhancement (TREE™) are advantageously used (for a discussion of 

RTVT\ see, International PCT application No. PCT/FR01/01366 published 
as International PCT application No. WO 01/186291 and the EXAMPLES 
below) or a refinement of that method provide herein and designated 
Tagged Replication and expression enhancement (TREE™) described 
5 above and in the examples. 
C. Practice of the process 

In one embodiment, the process provided herein includes the 
following steps. 

1 . Generation of diversity or source of existing diversity 

Generation of a plasmid library containing the genetic variants. The 
genetic variants are physically separated from each other. Any model 
such as, but not limited to, amino acid scanning, mutagenesis, or 
recombination may be used to generate the plasmid library. 

2. Expression of the genetic variants 
Any method for expression of variants is contemplated. In 

particular the following alternatives are particularly suitable for high 
throughput performance. 

a. Expression in bacterial hosts 
The mutated forms of the nucleic acid are prepared or introduced 
into plasmids for expression in bacterial cells. The genetic variants are 
expressed from suitable bacterial cells, which are prepared by 
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transformation aliquots of the cells with each member of the plasmid 
library (each genetic variant continues to be physically separated from 
each other). 

b. Expression in eukaryotic host cells 

5 A virus library is generated from the plasmid library. The virus 

library, in which each different member is separately maintained, is 
prepared by: 

(1) Transfection of the plasmid library into 
appropriate virus-producer cells (viruses produced, each one carrying a 

0 different genetic variant present in the original plasmid library, are 
physically separated from each other); 

(2) Titration of the virus library (of each individual 
virus present in the library, separately). Titration is effected by any 
method, but generally by either a method designated Real Time Virus 

5 Titering (RTVT™) (see, International PCT application No. 

PCT/FR01/01366 published as International PCT application No. 
WO 01/186291 and the EXAMPLES below) or a refinement of that 
method provide herein and designated Tagged Replication and expression 
enhancement (TREE™) described above and in the examples; 

(3) Standardization of the virus library to equal 
concentrations of all the individual viruses in the library (individual viruses 
continue to be physically separated from each other); 

(4) Expression of the genetic variants from 
appropriate mammalian cells by transduction with the virus library (each 
genetic variant continues to be physically separated from each other and 
each individual virus is handled separately from the others). 

3. Phenotypic characterization of the variant proteins. 
The variant proteins are expressed (from either plasmids in bacterial 
cells (step 2) or viruses in mammalian cells (step 4)) and their activity is 
assessed in one or more appropriate specific assays. The assays can be 
both types: biochemical (cell-free) assays and/or cell-based assays. The 
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variant proteins in the library are physically separated from each other and 
their activity is individually assessed on a one-by-one basis. 

The assays can be performed in one of a variety of ways, 
including, but are not limited to: 
5 a. Using a single-point dilution for each individual variant 

protein, followed by a kinetic analysis (multiple time points) of the read- 
out by technologies like Tagged Replication and expression enhancement 
(TREE™), or any other appropriate technology 

b. Using serial dilutions of each individual variant protein, 
5 10 followed by, for example, the Hill-based analysis of the read-out by 

technologies or any other appropriate technology. Hill based analyses 
i! assess the interaction between cells and biological agents (see, published 

International PCT application No. WO 01/44809 based on PCT n° 
PCT/FROO/03503). 

15 The goal of these methods is to identify proteins having an evolved 

function or property. 

Lead identification 
Based on the results obtained from the assays described above, 
each individual protein variant is individually tested for the parameters 
20 that assess the activity, property, function or structure of interest. 

Variants are ranked out according to their activity features. Those variant 
proteins best suited for the specificities of each individual project and 
system under study are then selected. The selected leads can be used for 
the desired purpose or further evolved or mutated to achieve desired 
25 activities. 

Typically, as for most directed evolution methods, the process is an 
iterative one, in which mutated variants are produced, screened, the best 
identified and then selected. The selected variants are then subjected to 
further evolution and the screening process repeated. This is repeated 
30 until the desired goal is achieved. 
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iher evolution may employ the methods herein or any 
lion method or combinations thereof. The methods for 
tion will include the amino-acid scan method herein, which 
onal approach to variant generation. Other rounds can 
lations of any other method for directed evolution known 
ations thereof, 
evolution of a viral gene 

inant viruses have been developed for use as gene therapy 

therapy applications are hampered by the need for 
f vectors with traits optimized for this application. The 
Jt methods provided herein are ideally suited for 
f such vectors. In addition to use for development of 
ral vectors for gene therapy, these methods can also be 
and modify the viral vector backbone architecture, trans- 
i helper functions, where appropriate, regulatable and 
promoters and transgene and genomic sequence analyses. 
,AV (rAAV) is a gene therapy vector that can serve as a 
cation of the methods herein for these and other purposes, 
ssociated virus (AAV) is a defective and non-pathogenic 

requires co-infection with either adenovirus or herpes 
ovide helper functions, for its growth and multiplication, 
ensive body of knowledge regarding AAV biology and 
?.g., Weitzman et al. (1996) J. Virol. 70: 2240-2248 

eta/. (1 997) J. Virol. 77:2722-2730; Urabeefa/. (1999) 
82-2693; Davis et al. (2000) J. Virol. 23:74:2936-2942; 
01) J. Virol. 75:3230-3239; Deng et al. (1992) Anal 
1-85; Drittanti et al. (2000) Gene Therapy 7:924-929; 
/. (1983) J. Virol. 45:555-564; Hermonat et al. (1984) J. 
:39; Chejanovsky et al. (1989) Virology 7 73:120-128; 

al. (1 990) J. Virol. 64:1764-1770; Owens et al. (1991) 
4-22; Owens et al. (1 992) J. Virol. 66:1236-1240; 
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Qicheng Yang et al. (1992) J. Virol. 66:6058-6069; Qicheng Yang et al. 
(1993) J. Virol. 67:4442-4447; Owens et al. (1993) J. Virol. 62:997- 
1005; Sirkka et al. (1994) J. Virol. 65:2947-2957; Ramesh et al. (1995) 
Biochem. Biophy. Res. Com. Vol 210 (3), 717-725; Sirkka (1995) J. 
Virol. 65:6787-6796; Sirkka et al. (1996) Biochem. Biophy. Res. Com. 
220:294-299; Ryan et al. (1996) J. l//>o/. 70:1542-1553; Weitzman et al. 
(1 996) J. Virol. 70:2440-2448; Walker et al. (1997)7. Virol. 77:2722- 
2730; Walker et al. (1997) J. Virol. 77:6996-7004; Davis ef a/. (1999) J. 
Virol. 73:2084-2093; Urabe et al. (1 999) J. Virol. 73:2682-2693; Gavin 
efaA (1999) J. Virol. 73:9433-9445; Davis et al. (2000) J. Virol. 
74:2936-2942; Pei Wu et al. (2000) J. WroA 74:8635-8647; Alessandro 
Marcello et al. (2000) J. Virol. 74:9090-9098). AAV are members of the 
family Parvoviridae and are assigned to the genus Dependovirus. 
Members of this genus are small, non-enveloped, icosahedral with linear 
and single-stranded DNA genomes, and have been isolated from many 
species ranging from insects to humans. 

AAV can either remain latent after integration into host chromatin 
or replicate following infection. Without co-infection, AAV can enter host 
cells and preferentially integrate at a specific site on the q arm of 
chromosome 19 in the human genome. 

The AAV genome contains 4975 nucleotides and the coding 
sequence is flanked by two inverted terminal repeats (ITRs) on either side 
that are the only sequences in cis required for viral assembly and 
replication. The ITRs contain palindromic sequences, which form a hairpin 
secondary structure, containing the viral origins of replication. The ITRs 
are organized in three segments: the Rep binding site (RBS), the terminal 
resolution site (TRS), and a spacer region separating the RBS from the 
TRS. 

Regulation of AAV genes is complex and involves positive and 
negative regulation of viral transcription. For example, the regulatory 
proteins Rep 78 and Rep 68 interact with viral promoters to establish a 
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feedback loop (Beaton et al. (1989) J. Virol 53:4450-4454; Hermonat 
(1994) Cancer Lett 37:129-136). Expression from the p5 and p19 
promoters is negatively regulated in trans by these proteins. Rep 78 and 
68, which are required for this regulation, have bind to inverted terminal 
repeats (ITRs; Ashktorab et al. (1989) J. Virol. 63:3034-3039) in a site- 
and stand-specific manner, in vivo and in vitro. This binding to ITRs 
induces a cleavage at the TRS and permits the replication of the hairpin 
structure, thus, illustrating the Rep helicase and endonuclease activities 
(Im et al. (1990) Cell 61:447-457; and Walker et al. (1997) J. Virol. 
77:6996-7004), and the role of these non-structural proteins in the initial 
steps of DNA replication (Hermonat et al. (1984) J. Virol. 52:329-339). 
Rep 52 and 40, the two minor forms of the Rep proteins, do not bind to 
ITRs and are dispensable for viral DNA replication and site-specific 
integration (Im et al. (1992) J. Virol. 66:1119-112834; Ni et al. (1994) J. 
Virol. 63:1128-1138. 

The genome (see, FIG. 4) is organized into two open reading 
frames (ORFs, designated left and right) that encode structural capsid 
proteins (Cap) and non-structural proteins (Rep). There are three 
promoters: p5 (from nucleotides 255 to 261: TATTTAA), p19 (from 
nucleotide 843 to 849: TATTTAA) and p40 (from nucleotides 1822 to 
1827: ATATAA). The right-side ORF (see FIG. 4) encodes three capsid 
structural proteins (Vp 1-3). These three proteins, which are encoded by 
overlapping DNA, result from differential splicing and the use of an 
unusual initiator codon (Cassinoti et al. (1988) Virology 767:176-184). 
Expression of the capsid genes is regulated by the p40 promoter. Capsid 
proteins VP1 , VP2 and VP3 intiate from the p40 promoter. VP1 uses an 
alternate splice acceptor at nucleotide 2201; whereas VP2 and VP3 are 
derived from the same transcription unit, but VP2 use an ACG triplet as 
an initiation codon upstream from the start of VP3. On the left side of 
the genome, two promoters p5 and p19 direct expression of four 
regulatory proteins. The left flanking sequence also uses a differential 
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splicing mechanism (Mendelson eta/. (1986) J. V7ro/ 60:823-832) to 
encode the Rep proteins, designated Rep 78, 68, 52 and 40 on the basis 
molecular weight. Rep 78 and 68 are translated from a transcript 
produced from the p5 promoter and are produced from the unspliced and 
5 spliced form, respectively, of the transcript. Rep 52 and 40 are the 
translation products of unspliced and spliced transcripts from the p19 
promoter. 

The rep protein is a adeno-associated virus protein involved in a 
number of biological processes necessary to AAV replication. The 
10 production of the rRep proteins enables viral DNA to replicate, 

encapsulate and integrate (McCarty eta/. (1992) J. Virol 66:4050-4057; 
Horer eta/. (1995) J. Virol 63:5485-5496, Berns eta/. (1996) Biology of 
Jh Adeno-associated virus, in Adeno-associated virus (AAV) Vectors in Gene 

Therapy, K.I. Berns and C. Giraud, Springer (1996); and Chiorini et al. 
fU 15 (1996) The Roles of AAV Rep Proteins in gene Expression and Targeted 

Integration, from Adeno-associated virus (AAV) Vectors in Gene Therapy, 
K.I. Berns and C. Giraud, Springer (1996)). A rep protein with improved 
activity could lead to increased amounts of virus progeny thus allowing 
higher productivity of rAAV vectors. 
20 AAV and rAAV have many applications, including use as a gene 

transfer vector, for introducing heterologous nucleic acid into cells and for 
genetic therapy. Advances in the production of high-titer rAAV stocks to 
the transition to human clinical trials have been made, but improvement of 
rAAV production will be complemented with special attention to clinical 
25 applications of rAAV vectors as successful gene therapy approach. 
Productivity of rAAV (i.e. the amount of vector particles that can be 
obtained per unitary manufacturing operation) is one of the rate limiting 
steps in the further development of rAAV as gene therapy vector. 
Methods for high throughput production and screening of rAAV have 
30 been developed (see, e.g., Drittanti et al. (2000) Gene Therapy 7:924- 
929) Briefly, as with the other steps in methods provided herein, the 
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plasmid preparation, transfection, virus productivity and titer and 
biological activity assessment are intended to be performed in 
automatable high throughput format, such as in a 96 well (or other 
number or multiples thereof, such as 384, 1536 . . . 9600, 9984 . . .) 
5 formats. 

Since the Rep protein is involved in replication it can serve as a 
target for increasing viral production. Since it has a variety of functions 
and its role in replication is complex, it has heretofore been difficult to 
identify mutations that result in increase viral production. The methods 
p 10 herein, which rely on in vivo screening methods, permit optimization of its 

activities as assessed by increases in viral production. Provided herein 
fU are Rep proteins and viruses and viral vectors containing the mutated Rep 

proteins that provide such increase. The amino acid positions on the rep 
proteins that are relevant for rep proteins activities in terms of AAV or 
15 rAAV virus production are provided. Those amino acid position are such 
that a change in the amino acid leads to a change in protein activity either 
to lower activity or increase activity. As shown herein, the alanine or 
amino acid scan revealed the amino acid positions important for such 
activity (i.e. hits). Subsequent mutations produced by systematically 
20 replacing the amino acids at the hit positions with the remaining 18 amino 
acids produced so-called "leads" that have amino acid changes and result 
in higher virus production. In this particular example, the method used 
included the following specific steps. 
Amino acid scan 

25 In order to first identify those amino acid (aa) positions on the rep 

protein that are involved in rep protein activity, an Ala-scan was 
performed on the rep sequence. For this, each aa in the rep protein 
sequence was individually changed to Alanine. Each resulting mutant rep 
protein was then expressed and the amount of virus it could produced 
30 measured as indicated below. The relative activity of each individual 
mutant compared to the native protein is indicated in FIG 3A. HITS are 
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those mutants that produce a decrease in the activity of the protein (in 
the example: all the mutants with activities below about 20 % of the 
native activity). 

In a second experimental round, which included a new set of 
mutations and phenotypic analysis, each amino acid position hit by the 
Ala-scan step, was mutated by amino acid replacement of the native 
amino acid by the remaining 18 amino acids, using site directed- 
mutagenesis. 

In both rounds, each mutant was individually designed, generated 
and processed separately, and optionally in parallel with the other 
mutants. Neither combinatorial generation of mutants nor mixtures 
thereof were used in any step of the method. 

A plasmid library was thus generated in which each plasmid contained 
a different mutant bearing a different amino acid at a different hit 
position. Again, each resulting mutant rep protein was then expressed 
and the amount of virus it could produced measure as indicated below. 
The relative activity of each individual mutant compared to the native 
protein is indicated in FIGURE 3B. LEADS are those mutants that lead to 
an increase in the activity of the protein (in the example: the ten mutants 
with activities higher, typically between 6 to 10 times, than the native 
activity). 

Expression of the genetic variants and phenotypic characterization. 

The rep protein acts as an intracellular protein through complex 
interaction with a molecular network composed by cellular proteins, DNA, 
AAV proteins and adenoviral proteins (note: some adenovirus proteins 
have to be present for the rep protein to work). The final outcome of the 
rep protein activity is the virus offspring composed by infectious rAAV 
particles. It can be expected that the activity of rep mutants would affect 
the titer of the rAAV virus coming out of the cells. 

As the phenotypic characterization of the rep variants can only be 
accomplished by assaying its activity from inside mammalian cells, a 
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mammalian cell-based expression system as well as a mammalian cell- 
based assay was used. The individual rep protein variants were expressed 
in human 293 HEK cells, by transfection of the individual plasmids 
constituting the diverse plasmid library. All necessary functions were 

provided as follows: 

(a) the cellular proteins present in the permissive specific 293 HEK 

cells; 

(b) the AAV necessary proteins and DNA were provided by co- 
transfection of the AAV cap gene as well as a rAAV plasmid vector 
providing the necessary signaling and substrate ITRs sequences; 

(c) the adenovirus (AV) proteins were provided by co-transfection 
with a plasmid expressing all the AV helper functions. 

A library of recombinant viruses with mutant rep encoding genes 
was generated. Each recombinant, upon introduction into a mammalian 
cell and expression resulted in production of rAAV infectious particles. 
The number of infectious particles produced by each recombinant was 
determined in order to assess the activity of the rep variant that had 
generated that amount of infectious particles. 

The number of infectious particles produced was determined in a 
cell-based assay in which the activity of a reporter gene, in the 
exemplified embodiment, the bacterial lacZ gene, or virus replication (Real 
time PCR) was performed to quantitatively assess the number of viruses. 
The limiting dilution (titer) for each virus preparation (each coming from a 
different rep variant) was determined by serial dilution of the viruses 
produced, followed by infection of appropriate cells (293 HEK or HeLa 
rep/cap 32 cells) with each dilution for each virus and then by 
measurement of the activity of the reporter gene for each dilution of each 
virus. Hill plots (NAUTSCAN 1 ") as described herein (published as 
International PCT application No. WO 01/44809 based on PCT n° 
PCT/FROO/03503, Dec, 2000; see EXAMPLES) or a second order 
polynomial function (Drittanti et al. (2000) Gene Ther. 7: 924-929) was 
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used to analyze the readout data and to calculate the virus titers. Briefly, 
the titer was calculated from the second order polynomial function by 
non-linear regression fitting of the experimental data. The point where 
the polynomial curve reaches its minimum is considered to be the titer of 
the rAAV preparation. A computer program for calculation of titers has 
been developed (see Drittanti et al. (2000) Gene Then 7: 924-929) to 
assess the minima. 

The TREE method described herein can be used to analyze the 
readout data and to calculate the virus titers. The results are shown in the 



n 10 EXAMPLE below. 



Comparison between results of full-length Hit position analysis 
reporter here and the literature 

The experiments identified a number of heretofore unknown 
mutation loci, which include the hits at positions: 4, 20, 22, 28, 32, 38, 

15 39, 54, 59, 124, 125, 127, 132, 140, 161, 163, 193, 196, 197, 221, 
228, 231, 234, 258, 260, 263, 264, 334, 335, 341, 342, 347, 350, 
354, 363, 364, 367, 370, 376, 381, 389, 407, 411, 414, 420, 421, 
422, 428, 429, 438, 440, 451, 460, 462, 484, 488, 495, 497, 498, 
499, 503, 51 1, 512, 516, 517 and 518 with reference to the amino 

20 acids in Rep78 and Rep 68. Rep 78 is encoded by nucleotides 321- 

2,186; Rep 68 is encoded by nucleotides 321-1906 and 2228-2252; Rep 
52 is encoded by nucleotides 993-2186, and Rep 40 is encoded by amino 
acids 993-1906 and 2228-2252 of wildtype AAV. 

Also among these are mutations that may have multiple effects. 

25 Since the Rep coding region is quite complex, some of the mutations may 
have several effects. Amino acids 542, 598, 600 and 601, which are in 
the Rep 68 and 40 intron region, are also in the coding region of Rep 78 
and 52. Codon 630 is in the coding region of Rep 68 and 40 and non 
coding region of Rep 78 and 52. 

30 Mutations at 10, 86, 101, 334 and 519 have been previously 

identified, and mutations, at loci 64, 74, 88, 175, 237, 250 and 429, but 
with different amino acid substitutions, have been previously reported. In 
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all instances, however, the known mutations reportedly decrease the 
activity of Rep proteins. Among mutations described herein, are 
mutations that result in increases in the activity the Rep function as 
assessed by detecting increased AAV production. 
5 Lead identification. 

Based on the results obtained from the assays described herein (i.e. 
titer of virus produced by each rep variant), each individual rep variant 
was assigned a specific activity. Those variant proteins displaying the 
highest titers were selected as leads and are used to produce rAAV. 

10 In further steps, rAAV and Rep proteins that contain a plurality of 

mutations based on the hits (see Table in the EXAMPLES, listing the hits 
and lead sites), are produced to produce rAAV and Rep proteins that 
have activity that is further optimized. Examples of such proteins and 
AAV containing such proteins are described in the EXAMPLES. 

15 The rAAV rep mutants are used as expression vectors, which, for 

example, can be used transiently for the production of recombinant AAV 
stocks. Alternatively, the recombinant plasmids may be used to generate 
stable packaging cell lines. To create a stable producer cell line, the 
recombinant vectors expressing the AAV with mutant rep genes, for 

20 example, are cotransfected into host cells with a plasmid expressing the 
neomycin phosphotransferase gene (neor) by transfection methods well 
known to those skilled in the art, followed by selection for G418 
resistance. 

Also among the uses of rAAV, particularly the high titer stocks 
25 produced herein, is gene therapy for the purpose of transferring genetic 
information into appropriate host cells for the management and correction 
of human diseases including inherited and acquired disorders such as 
cancer and AIDS. The rAAV can be administered to a patient at 
therapeutically effective doses. A therapeutically effective dose refers to 
30 that amount of the compound sufficient to result in amelioration of 
symptoms of disease. 
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Gene therapy 

Toxicity and therapeutic efficacy of the rAAV can be determined by 
standard pharmaceutical procedures in cell cultures or experimental 
animals, e.g., for determining the LDS 50 (the dose lethal to 50% of the 
5 population) and the ED50 (the dose therapeutically effective in 50% of 
the population). The dose ratio between toxic and therapeutic effects is 
the therapeutic index and it can be expressed as the ratio LD 50 /ED 50 . 
Doses that exhibit large therapeutic indices are preferred. Doses that 
exhibit toxic side effects may be used, care should be taken to design a 

10 delivery system that targets rAAV to the site of treatment in order to 
minimize damage to untreated cells and reduce side effects. 

The data obtained from cell culture assays and animal studies can 
be used in formulating a range of dosage for use in humans. The dosage 
of such rAAV lies preferably within a range of circulating concentrations 

15 that include the ED 50 with little or no toxicity. The dosage may vary 
within this range depending upon the dosage form employed and the 
route of administration utilized. A therapeutically effective dose can be 
estimated initially from cell culture assays. A dose may be formulated in 
animal models to achieve a circulating plasma concentration range that 

20 includes the IC 50 (ie., the concentration of the test compound which 
achieves a half-maximal infection or a half- maximal inhibition) as 
determined in cell culture. Such information can be used to more 
accurately determine useful doses in humans. Levels in plasma may be 
measured, for example, by high performance liquid chromatography. 

25 

The following examples are included for illustrative purposes only 
and are not intended to limit the scope of the invention. The specific 



-55- 



methods exemplified can be practiced with other species. The examples 
are intended to exemplify generic processes. 

EXAMPLE 1 

Titering or assessment of concentration by a method designated Real 
Time Vector Titering (RTVT™) 

This Example is based on the method described in International PCT 
application No. PCT/FR01/01 366, based on French application 
n° 0005852, filed 9 May 2000, and published as International PCT 
application No. WO 01/186291. This method assesses the 
titer or concentration of a biological agent (virus, gene transfer vector) in 
a sample, by measurement of the kinetics of change of a reporter 
parameter following the exposure of cells to the biological agent. 

As noted above, reporter parameters may include, but are not 
limited to,: gene / transgene expression related to the gene/transgene 
products, such as enzymatic activity, fluorescence, luminescence, antigen 
activity, binding to receptors or antibodies, and regulation of gene 
expression), differential gene expression, viral/vector progeny 
productivity, toxicity, cytotoxicity, cell proliferation and/or differentiation 
activity, anti-viral activity, morphogenetic activity, pathogenetic activity, 
therapeutic activity, tumor suppressor activity, oncogenetic activity, 
pharmacological activity. 

Serial dilution methods 

The assessment of the concentration or titer of biological agents 
using current approaches needs for serial dilutions of the agent. 
Serial dilutions of the agent are applied to a cell-based reported system, 
that elicits an output signal in response to the exposure to the agent. The 
intensity of the signal is a function of the concentration of the agent. The 
titer or concentration of the agent is determined as the highest dilution 
that still elicits a measurable response in the output. The higher the 
number of dilutions tested, the higher the accuracy of the value obtained 
for the titer. 
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This approach requires a set of serial dilutions for every biological 

agent whose titer needs to be determined. Thus, the application of this 

approach to the simultaneous titration of a large number of different 

biological agents is limited by the number of experimental points needed 

(example: for 30 biological samples: 20 serial dilutions x 30 biological 

agents: 600 experimental points). 

The approach in International PCT application No. 
PCT/FR0 1/0 1366 published as International PCT application 
No. WO 01/186291 

The intensity of the output signal (after exposure of reporter cells 
to the biological agent) is not only dependent on the concentration of the 
agent but also on the time after exposure. As time increases, the intensity 
of the signal increases. The kinetics of change of the intensity over time 
depends upon the concentration of the agent. Thus, lower concentrations 
of the agent will require longer times for the intensity to reach a given 
value that would be reached in shorter times after exposure to higher 
concentrations of the same agent. 

This approach (designated Real Time Virus Titering (RTVT"*) uses 
the following: a reference plot representing the relationship between the 
concentration of the agent and the time necessary for the intensity to 
reach a given threshold value is obtained using a reference preparation of 
biological agent, whose concentration or titer is known. This plot is then 
used to obtain the concentration of the biological agent under study by 
entering the time that a dilution of that agent needed for the intensity to 
reach the threshold value. 

Using this approach, there is no need for serial dilutions of the 
biological agent(s) under study. Once the reference plot (t/? vs c) is 
obtained, it can be used for the determination of the concentration or titer 
of any number of biological agents. Only one dilution of the biological 
agent under study is necessary to obtain the corresponding value of t{3 
that is then used to obtain the concentration or titer using the reference 
plot. 
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Thus, the application of this approach to the simultaneous titration 
of large numbers of different biological agents is facilitated by the fact 
that only one dilution of each sample is needed (example: for 30 
biological samples: 1 dilution x 30 biological agents: 30 experimental 
points (compared to 600 needed with the current approach). 
This approach is specially suited for the high throughput assessment of 
concentration or titer of large numbers of biological agents. 

The system 

The system includes the following elements: 

a preparation of the biological agent (virus, gene transfer vector, 
protein,...) whose concentration or titer is unknown and has to be 
determined. 

a reporter system including culture of a cell line (or a mixture 
of cell lines) that reacts to the exposure to the biological 
agent by displaying a specific output signal, 
a master preparation of a reference biological agent, of 
known concentration or titer, that is able to generate the 
output signal when the reporter cells are exposed to it. 
Practice of the method 

When the reporter cells are exposed to either the biological agent 
under study or the reference biological agent, an output signal is 
generated, that can be measured. 

The intensity of the output signal is called /; the concentration of the 
biological agent used is called c, the time of exposure of the cells to the 
biological agent is called t. The intensity of the output signal (/ ) is a 
function of c and t : 

/ increases as the concentration (c ) of the biological agent applied 

to the cells increases; 

/ increases as the time of exposure of the cells to the biological 
agent it) increases. 
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If the time t after exposure of the cells to the biological agent is 
kept constant, then, / will change as a direct function of c . 

If the concentration c of the biological agent is kept constant, 
then, / will be a direct function of f. 

fi is defined here as a threshold value of the intensity of the output 
signal, arbitrarily defined for every system under study. 

tfi is defined here as the time necessary to reach the threshold fi. 

Use of fi and tfi to determine the concentration or titer 
of a biological agent. 

The reporter cells are exposed to serial dilutions of a reference 
biological agent, whose concentration (or titer) is previously known. The 
intensity of the output signal (#) is measured at several time points (t) for 
every concentration (serial dilutions) c of the reference biological agent. 

/ is plotted vs t , and that, for every concentration c used of the 
reference biological agent. 

Using the plots obtained above, and for every concentration c of 
the reference biological agent, the time (tfi) necessary for the intensity of 
the output signal to reach a threshold value fi is obtained. 

With the data obtained above, tfi is plotted vs c. 

This plot represents the time necessary for the intensity of the 
output signal to reach the threshold value fi as a function of the 
concentration of the biological agent used. T his is a standard plot and 
will be used to determine the unknown concentration of the biological 
agent under study by measuring the time that a given dilution of it needs 
to give an output signal whose intensity equals the threshold fi. 

The reporter cells are exposed to a dilution of the biological agent 
under study (whose concentration or titer is to be determined). The 
intensity of the output signal (/) is measured over time until it reaches the 
threshold value fi. The time necessary for / to reach the value fi is 
recorded as tfi. 

The tfi value recorded above is entered into the standard plot 
obtained above) and the corresponding c value is obtained. 

-59- 



This c value represents the concentration or titer of the biological 
agent under study. 

Example of the Real Time Virus Titerng RTVT" 
titering method 

Rat-2 cells were infected with serial dilutions of a reference 

preparation of a retroviral vector carrying the green fluorescent protein 

(GFP) gene (vector pSI-EGFPI see, Ropp et al. (1995) Cytometry 27:309- 

317). At increasing times after infection, the level of expression of the 

transgene was determined (as the level of fluorescence due to the GFP 



gene) as the output signal. 

Table 3 represents the values obtained: 



Concentration 


Time after infection 


Output signal 


(1 = 10 6 particles/ml) 


(hrs) 


fluorescence 


0.1 


16 


20.4 




24 


30.1 




40 


95.1 




48 


138.7 




64 


157.3 


0.25 


16 


26.8 




24 


48.5 




40 


173.3 




48 


228.2 




64 


191.7 


0.5 


16 


38.1 




24 


72 




40 


198.7 




48 


296.2 




64 


203.7 



The threshold value of /? = 100 was arbitrarily selected for this example. 
The time (t/J) necessary for the output signal to reach the threshold /? , for 
every concentration is shown in table 4. 



Table 4 



Concentration 
(1 = 10 6 particles/ml) 


(hrs) 


0.1 


42 


0.25 


31 
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A plot of tfi versus concentration for the reference virus shows that 
the concentration and tfi exhibit a clearly defined relationship, that allows 
for the calculation of the concentration (c) of a sample, if the 
corresponding tfi of that sample is known. 

EXAMPLE 2 

* 

Tagged Replication and Expression Enhancement (TREE) for titering 

As discussed above, TREE is a method for titering and 
standardization of preparations of viruses, vectors, antibodies, libraries, 
proteins, genes and any other moiety that is detectable based upon a 
output signal, such as fluorescence. The TREE method is an improvement 
of the Real Time Virus Titering (RTVT) method (see, International PCT 
application No. PCT/FR01/01 366 published as International PCT 
application No. WO 01/186291). It is performed a reporter moiety, such 
as a reporter virus (with a known titer) and the test sample (with 
unknown titer). The reporter, such as a reporter virus has a readily 
detectable output signal that can be measured as a function of time. The 
effect of the moiety, such as a virus, of unknown titer is assessed. The 
moiety whose titer is assessed either increases or decreases the output 
signal as a function of time. This change in signal is used to assess the 
amount or concentration of the moiety of unknown concentration, and 
hence its titer. 

The method is exemplified herein using an AAV system for the 
determination of the titer of an AAV vector and an AAV-reporter vector 
as a competitor and wild type Adenovirus as helper virus. One of skill in 
the art readily can adapt the method to other systems, including other 
viruses, and other moieties for which a reporter system can be developed. 
Other such moieties include, but are not limited to, viral vectors, 
plasmids, libraries, proteins, antibodies, vaccines, genes, and nucleic acid 
molecules. 
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Materials and Methods 

1 . Cells and Viruses 

HeLa rep-cap32 cells, a HeLa derived cell line (kindly provided by P. 
Moullier, Laboratoire de Therapie Genique, CHU, Nantes; see, Salvetti et 
5 aL (1998) Hum Gene Ther 20:695-706; Chadeuf et al. (2000) J Gene 
Med 2:260-268) was grown in DMEM with 10% fetal calf serum. These 
cells were plated 24 h before infection at a density of 1 x 10 4 cells in 
single well of 96-well plates. rAAV-LacZ (10 10 ip/ml), rAAV-GFP (10 9 
ip/ml) vectors and Human Adenovirus type 5 (Ad5) (10 11 pfu/ml) were 

10 from CHU, Nantes. 

Hela rep-cap32 cells had been produced by cotransfecting plasmid 
pspRC, which harbors the AAV rep-cap genome with the ITRs deleted (bp 
190 to 4484 of wild-type AAV), with plasmid PGK-Neo, conferring 
resistance to G418 on Hela cells (see, Chadeuf et al. (2000) J. Gene 

15 Med. 2:260-268 and Salvetti et aL (1998( Hum Gene Ther. 3:695-706). 
Hela rep-cap 32 cells are a packaging line that harbor one copy of the 
genome with the ITRs deleted (see, also Tessier et aL (2001) J. Virol. 
75:375-383). 

Plasmid pspRC contains the AAV genome (positions 190-4,484 bp) 
20 with the ITRs deleted and was obtained by excising the rep-cap fragment 
(Xba\ fragment) from the well-known vector psub201 (Samulski et al. 
(1987) J Virol 67:3096-3101; also called pSSV9) by Xba\ digestion and 
inserting it into the Xba\ site of plasmid pSP72 (Promega). Plasmid 
psub201 (see, e.g., U.S. Patent No. 5,753,500) is a modified full-length 
25 AAV type 2 genomic clone contains all of the AAV type 2 wild-type 
coding regions and cis acting terminal repeats. 

2. Infection and measurement 

Four serial dilutions of a rAAV-LacZ (0.01, 0.0075, 0.005 and 
0.0025 see Table 2 below, designated samples 1-4, respectively) were 
30 made and used for co-infection of HeLa rep-cap32 cells together with 8 
different Ad5 multiplicity of infection (MOI; from 0.1 to 100/cell) and with 
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10 3 ml (10 6 infectious particles (ip)) or 10 4 ml (10 5 ip) rAAV-GFP viral 
vector. All the samples were done in triplicate. After infection, the plates 
were read at different times, from 34.5 h to 80 h (every 30 minutes). 

rAAV-GFP is an SSV9-derived vector; SSV9 is a clone containing 
the entire adeno-associated virus (AAV) genome inserted into the Pvull 
site of plasmid pEMBL (see, Du eta/. (1996) Gene Ther 3:254-261 ). The 
rAAV-GFP and rAAV-LacZ plasmids are SSV9 with a GFP or LacZ gene 
under control of the cytomegalovirus (CMV) immediate-early promoter. 
All the samples were done in triplicate. 

3- Process 

Figure 2 shows the overall procedure in 96 well format. Cells were 
plated 24 h before infection. Co-infection of rAAV-GFP with serial 
dilutions of rAAV-LacZ together with Ad5 (different MOD, were done. 
Then the plates were read at different times using the Analyst AD&HT 
micro plate-reader (LJL BioSystems). 

4. Analysis 

For this kinetic technique, Fluorescence Intensity (Fl) of the infected 
cells is measured as a function of the time. Serial dilutions of the AAV- 
competitor vector AAV-lacZ vector, which decreases the fluorescence 
signal, are performed. For this example, fluorescence was measured for 
AAV-GFP with 10 6 ip and 10 5 ip and then 10 6 ip of the AAV-GFP reporter 
virus in the serial dilutions of the competitor virus, AAV-lacZ vector in a 
96-well format (samples 1-4, see Table below). Measurements were 
taken of each well and curves of Fl (of the GFP) versus time (hrs) were 
obtained (see Figure 2B). 

An arbitrary one value for Fl (see Fig. 2B, 6 x 10 6 Fl units), 
typically, though not necessarily, near the greatest separation among the 
curves so that the numbers are readily discernable, was selected. The 
point at which each of the curves intersect this value is beta time (t/?) for 
that combination of amounts of reporter plus dilution of the virus of 
unknown titer, tjff, taken from the Fl vs. Time (hrs) curves, for each 
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sample containing a dilution of the unknown plus 10 6 ip of the reporter 
virus is set forth in column 2 of Table 2 below. 

To determine the titer of the test virus, the t£ for the AAV GFP 
(reporter virus) is plotted versus quantity of ip (i.e a straight line between 
the tjff for the 10 6 ip and the 10 5 ip) (Fig. 2C). For any tfi of the unknown 
virus, the quantity of ip can be determined from this curve. The beta 
time {t/3) of each sample (in this case for the different dilutions of rAAV 
LacZ mixed with 10 6 infectious particles of rAAV-GFP) is determined, and 
then the residual number of infectious particles of rAAV-GFP for each 
sample. The difference between 10 6 ip of rAAV-GFP put in each sample 
and the number of ip detected by fluorescence in the same well is the 
actual quantity of rAAV-GFP competed (consumed) by the unknown 
rAAV (in this case rAAV-LacZ). This number is determined for each 
dilution. The quantity rAAV-GFP consumed is the same quantity of 
unknown rAAV in the sample. This quantity is present in one volume of 
unknown rAAV, which in this example is 1 ml. Based upon this, the 
infectious titer of the unknown rAAV is determined. The results are 
shown in Table 2. 

TABLE 2 

AAV LacZ titration by TREE titration 



Sample 


volume 


Xfi (hrs) 


Residual 
AAV-GFP 


Consumed 
AAV-GFP 


AAV LacZ 
Concentration 
(i.p./10 2 /vl) 


1 


10.0 x 10 3 


66.5 


5.56 x 10 5 


4.44 x 10 5 


4.44 x 10 5 


2 


7.5 x 10 3 


66.5 


5.56 x 10 5 


4.44 x 10 5 


5.93 x 10 5 


3 


5.0 x 10' 3 


65 


7.06 x 10 5 


2.94 x 10 5 


5.88 x 10 5 


4 


2.5 x 10* 3 


64 


8.06 x 10 5 


1.94 x 10 5 


7.76 x 10 5 



The average titer using this method was 6.01 x 10 ,u ip/ml (6.0 



x 10 7 



ip//vl = 6.01 x 10 5 ip/0.01 The standard deviation was 1.37 x 10 



10 



ip/ml with an error of ±23%. 
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EXAMPLE 3 

Hill Analysis of the sere ning assay output 

It is important to have reliable methods for screening and/or 
evaluating the performance of a set of biological agents, such as a library 
5 of viral or non-viral recombinant vectors, vaccines, recombinant proteins 
and antibodies, in a complex biological system, such as living target cells 
When developing such agents, for example gene therapy vectors and 
other agents for therapeutic use, it is necessary to be able to evaluate 
and compare performance among candidates. 
10 The progress of gene transfer into gene therapy depends upon the 

rij capacity to develop gene transfer vectors into therapeutic drugs. Clinically 

Til 

fij relevant vectors need to be efficient and safe, in reaching and infecting 

target cells and in ensuring a persistent level of expression of the 
therapeutic gene with a minimum of adverse effects. The availability of 
15 standardized quantitative methods, suitable for an accurate and objective 
assessment of titer, performance and safety, is necessary for the 
O pharmaceutical development of gene vectors as drugs. 

Any method for assessment is contemplated herein as long as it is 
adapted for use in a high throughput format. Of particular interest is the 
20 Hill equation based method of International PCT application No. 

WO 01/44809 (International PCT application No. PCT/FR00/03503, based 
on French application FR 9915884). 

Two widely used parameters that provide quantitative information 
about the potential performance of a gene transfer vector preparation are 
25 the titer of physical particles and the titer of infectious particles. Vector 
preparations with high titer of infectious particles and low physical 
particles/infectious particles ratio are considered to be of higher quality. 

The titer in physical particles(pp) (see, e.g., Mittereder eta/. (1996) 
J, Virol. 70:7498-7509; Atkinson et al. (1998) Nucl. Acids Res. 26:2821- 
30 2823; Kechli et al. (1998) Hum. Gene Ther. 5:587-590; and Nelson et al. 
(1998) Hum. Gene Ther. 3:2401-2405), which represents the total 
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number of vector particles, is usually evaluated from the vector content 
by detecting the nucleic acid contents (nucleic acids hybridization and 
OD 260 respectively for AAV and AdV), detecting viral protein content (for 
example, reverse transcriptase (RT) activity and p24 content for MLV and 
5 HIV, respectively). 

Among the physical particles (pp), there are particles potentially 
active in performing transduction (ip, infectious particles), as well as 
particles that are inactive (nip, non-infectious particles) (Ruffing et al. 
(1994) J.Gen.ViroL 75:3385-3392; Kechli et al. (1998) Hum. Gene Ther. 

10 5:587-590.). The pp and the ip/nip ratio, are features of the packaging 
system, the manufacturing process and the vector itself. 

The infectious particles (ip) (infectious units, transducing units, 
etc.) are evidenced by the changes observed in the infected cells (vector 
DNA replication, provirus integration, cell lysis, transgene expression and 

15 other observable parameters. Infectious particles (ip) measures the 
number of particles effective in performing a process whose output is 
being measured; not all particles participate or are capable of participating 
in all processes. 

The precise assessment of ip is not straightforward. Existing 

20 methods are mainly based on serial dilution experiments followed by 
either linear extrapolation or asymptote approximates. The titer of 
infectious particles (ip: infectious unity, transduction unity) (see, e.g., 
Mittereder et al. (1 996) J. Virol. 70:7498-7509; Weitzman et al. (1996) 
J. Virol. 70: 1845-1 854; Salvetti et al. (1998) Hum. Gene Ther. 5:695- 

25 706) is evaluated by the studying observed changes in infected cells, 

such as viral replication, provirus integration, cellular lysis and transgene 
expression, using methods based on serial dilutions, followed either by a 
linear extrapolation or an asymptotic approximation. Thus, ip measures 
the number of active particles in the measured process; it includes 

30 physical particle (pp) and inactive particles (nip or non-infectious 
particles) . 
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In order to resolve the problem of the titer determination and the 
comparison of different recombinant viruses used in gene therapy, the 
variation of the particles/infectious power ratio has been used (see, e.g., 
Atkinson et at. (1998) Nucl. Acids Res. 25:2821-2823; and International 
PCT application No. WO 99/1 1764, which describe a method that uses 
step of amplification viral genetic material in a host cellular line, 
preparation of vectors of unknown titer obtained by serial dilution and an 
internal check of known titer). In particular, the method uses cells 
infected with a viral preparation in the different wells of a microtiter plate, 
viral genome replication in the host cells, nucleic acid hybridization and 
determination of the relative amount of replicated viral nucleic acid in 
each well. 

All of these methods measure the physical particles (pp) titer and/or 
measure infectious particles Up) titer in order to evaluate a gene transfer 
vector. A high quality vector preparation is one with an high titer of 
infectious particles and a low pplip ratio. These parameters provide 
quantitative information on the performance of a gene transfer vector. 
Because of the inaccuracy of the procedures used for assessing pp and 
especially ip, these parameters are not informative enough to precisely 
define the features of a gene therapy vector nor those of a particular 
preparation thereof. The actual procedures used for pp and ip evaluation 
change with the vector type, are not very reproducible nor exact, so 
these parameters do not contain enough information to allow a very fine 
definition of vector characteristics and performances. 

Hill equation-based analyses 

In this method complex biological processes, including those 
involving the response of cells (in vitro and in vivo) to biological agents, 
such as, for example, cells, viruses, vaccines, viral and non-viral gene 
vectors, antibodies, antigens, proteins in general and plasmids, are 
characterized using the formal analysis first introduced by Hill (see, Hill 
(1910) J. Physiol. 40AP; Hill (1913) I. Biochem.J. 7:471-480). 
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International PCT application No. WO 01/44809 (based on 
PCT/FR00/03503, priority claimed to French application FR 9915884) 
describes the use of the Hill equation (see, Hill (1910) J. Physiol. 40:4P; 
Hill (1913) /. Biochem.J. 7:471-480; see, International PCT application 
No. PCT/FR00/03503) for analysis and characterization of the biological 
and/or pharmacological activity of biological agents (viruses, vectors or 
cells) on biological assay systems in vitro (cell-based) or in vivo. 

A number of useful parameters, derived from the Hill equation, are 
scored and used to quantify relevant features of the biological agent, of 
the cells, as well as of the biological process or reaction involved. 

In particular, methods for calculation and analysis of the 
parameters of biological and pharmacological activity of native, 
attenuated, modified or recombinant viruses, vaccines, recombinant viral 
and non-viral gene transfer vectors, cells, antibodies and protein factors 
in in vitro (cell-based) or in vivo assays are described. This method is 
adapted for high throughput processes and is sufficiently accurate to 
allow a very fine definition of vector characteristics and performances. 

International PCT application No. WO 01/44809 provides, a 
standard process for evaluating the interaction between any biological 
agent, such as a gene therapy vector, with a complex biological system 
(living target cells). It provides a screening process for a pool of complex 
biological agents, in order to select test agents that have a desired 
property, activity, structure or whatever is being sought. 

Different biological agents and assay systems (cells) are compared 
and ranked out on the basis of their performance, assessed through the 
Hill parameters. Thus, the accurate analysis and comparison of the 
biological response of complex assay systems (in vitro and in vivo) to 
complex biological agents is achieved experimentally. The Hill-based 
analysis (rT,K,T,€,rj t 0) is used for a variety of purposes, including, but not 
limited to: 
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i) validation and optimization of the manufacturing processes used 
to obtain the biological agents; 

ii) development and optimization of the components of the 
biological agents (proteins, genomes, genetic units); 

iii) development and optimization of assays and analytical tests for 
the characterization of the biological agents. 

The method includes the steps of: 

(a) preparation, for each biological agent, of a sample scale, 
obtained by a serial dilution of the biological agent at a R1 concentration, 

(b) incubation of each sample of the dilution scale obtained in 1, 
with the target cells at a constant concentration R2, 

(c) determination of the P product from the reaction R1 + R2, at a 
t moment, in each the sample; and 

(d) realization of a theoretical curve H from the experimental points 
R1 and P, for each biological agent by iterative approximation of 
parameters of the reaction R1 + R2 P, at the t moment, in accordance 
with this equation: 

P = p max (A7R1) r / (k + (/7R1) r ) r=1 n (2) 

in which: 

R1 represents the biological agent concentration in a sample 
from the scale; 

R2 is concentration of target cells (in vitro or in vivo) 

P (output) represents the product from the reaction R1 + 
R2 at a t moment; 

P max represents the reaction maximal capacity; 

k represents, at a constant R2 concentration, the resistance 
of the biological system for responding to the biological agent (resistance 
constant R2); 

r represents a coefficient that depends on R1 and 
corresponds to the Hill coefficient; and 
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/7 represents the intrinsic power of the R1 biological agent to 
induce a response in the biological system (P production at the t 

moment), and 

(e) sorting the k and n values obtained in (d) for each 
5 biological agent and the biological agent, and then ranking according to 
the values thereof. 

Using the parameters {n,K,T,0 f €,n) the activity of a biological agent 
on a complex biological system, as well as its intrinsic features can be 
fully characterized and compared. In addition, different biological 
10 systems either in vitro cell-based) or in vivo can be compared. 

Hill Equation 
The Hill equation: 



■F\ r = n 

P = I P max . R1 r / (K 4- (R1) r ) R2 constant (1) 

y. 15 r = 1 



where R1, P, P max and K represent, respectively, the concentration 
of the reagent R1, the concentration of the product, the maximal capacity 
p, of the reaction and the 'affinity' constant between R1 and R2. The Hill's 

coefficient (r) is a function of R1 . The coefficient r is equal to 1 when 
20 independent non-interactive binding sites are involved between R1 and 
R2, such as in reactions that follow kinetics described by Michaelis- 
Menten; and r varies from 1 to n for systems where the sites involved in 
the interaction between the R1 and R2 are not independent from each 
other, and the affinity for R1 at any R2 binding site varies as a function of 
25 either i) the degree of occupancy of other R2 sites; ii) the concentration 
of R1 itself or iii) the concentration of other (positive or negative) 
regulators. 

The Hill equation, thus, is a general formalization that describes the 
interaction reaction between molecules. It expresses the amount of 
30 product formed as a function of the concentration of the reagents and of 
the affinity constant of the system. Originally developed for the study of 
the dissociation between haemoglobin and oxygen, the Hill equation 
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covers the formal Michaelis-Menten analysis of enzyme kinetics, the 

analysis of ligand-receptor binding and of the allosteric protein systems. 

According to Hill, for a simple reaction like 

R1 + R2 P 
K 

where the affinity between R1 and R2 changes with concentration 
of each, the Hill equation describes the accumulation of the product P as 
a function of the concentration of one of the reagents (R1) and of the 
intrinsic properties (K) of the system. 

This equation can be applied to complex biological systems. For 
example, the response of the cells to infection (P), can be analyzed by 
applying an Hill-type equation. The amount of cells growing in vitro (R2) 
are infected with increasing concentrations of recombinant viruses (R1), 
and (P) is monitored. A Hill equation is iteratively fitted to the 
experimental data. 

For analyses of viral output as exemplified herein, 

virus + cell -> transduced cell -> output (viral genome replication), 
Equation (1) is specifically reformulated as: 

P = p max (/7R1) r / (/c + (77R1) r ) r=1,...,n (2) 

where P, P max , R1, n, r and k , as described above, represent, 
respectively, the output signal (P) (the level of viral gene expression, or 
the level of virus replication), the maximal output signal (Pmax), the initial 
concentration (R1) of infectious viral particles (those susceptible to trigger 
the process leading to P), the potency of the vector (n; a factor that 
affects the concentration of the vector (R1) by its specific strength or 
activity, the Hill's coefficient (r) and the constant of resistance of the 
reaction or process (k). 
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The constant of resistance k 

The concept of k is analogous to those of dissociation, kinetics, 
equilibrium or affinity constants concepts for simple chemical and 
biological reactions, k is a feature of the process (reaction) and of the 
biological system tested (cell type), k is a key parameter for the 
characterization of the assay system and the assessment of its 
performance as a test for the reaction under study. 

k measures the internal resistance offered by the process or 
reaction triggered by the biological agent, to proceed to P. k is specific to 
a particular process or reaction tested. In addition k is specific to the 
particular biological system tested. Different cell lines and types will 
display different k for the same reaction. Moreover, factors affecting the 
performance of a cell to accomplish the reaction (like contaminants, toxic 
agents, etc.) affect k in that cell. 

Variations in k affect equation (2) by shifting the curve to the right 
or to the left, according to whether the value of k increases or decreases, 
respectively. All curves differing only in k are parallel each other. 
k finds its direct and practical application in i) assay development and 
validation and ii) assessment of the susceptibility or sensitivity of different 
cell types or tissues to undertake the reaction under study and to be 
affected by it. 

The potency n 

n measures the intrinsic potency of the biological agent to 
accomplish P against the resistance (k) offered by the reaction process. 
For every infectious virus particle (Rl) added to the assay, the actual 
activity of the virus added is given by /7R1 . In order to report an output P, 
the potency n has to push forward the reaction inside the cell against k. 
n is specific to the particular biological agent for the reaction under study. 
n is a feature of the biological agent. 
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Different versions or variants of the biological agent will display 
different rr for the same reaction. Thus, mutations, conformational 
changes or other modifications on the biological agent are expected to 
change its rr for a given reaction process. 

The concept of rr is analogous to that of chemical activity by 
opposite to concentration for simple compounds. 77 is a correction factor 
that affects the concentration (R1) of the biological agent to indicate its 
actual strength or activity for a given reaction process. 

Variations in rr affect equation (2) by shifting the curve to the right 
or to the left, according to whether the value of rr decreases or increases, 
respectively. Curves differing only in rr are not parallel each other. The 
slope of the curve given by equation (2) increases as rr increases. 

rr is a key parameter for the characterization of the biological agent 
and the assessment of its performance to accomplish the reaction under 
study. /7 finds its direct and practical application in i) biological agent 
optimization and development as it allows to compare the relative potency 
of variants of the agent. 

/7 is a valuable tool in the field of vaccine, gene transfer vector and 
antibody development, for the comparison between two or more different 
agents or different versions of the same agent, for performance. Two 
agents, for instance, may elicit equivalent potencies for gene transfer, 
while their potencies for immunogenicity be different. The use of rr t a 
quantitative and accurate parameter for assessing potency, will allow for 
ranking the candidates according to their potency (i.e., for gene transfer, 
gene expression, immunogenicity and other such properties and activities) 
and to make rational decisions about the relative value of the agent leads. 

The efficiency € 

€ measures the maximal global efficiency of the reaction process 
when a biological agent characterized by a given rr value interacts with a 
biological system characterized by a given k, e is specific to the particular 
couple biological agent (rr) I biological system (k) for the reaction under 
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study, e is a feature of the global reaction process and intervening 
reagents. Changes in either n , k, or both, will lead to changes in e. 

The efficiency of the reaction process described by equation (2) is 
given by the increase in the output P that can be obtained by increasing 
the input R1 . Thus, the first derivative of P with respect to R1, or the 
slope of the curve described by equation (2), gives the global efficiency of 
the reaction at every R1 input. The maximal global efficiency, or e, is 
given directly by either the slope at the inflection point of the curve 
described by equation (2) or by the maximum of its derivative <5P/<JR1. 
The slope of the curve given by equation (2) and the maximum of <$P/£R1 
increase as e increases. 

€ is a key parameter for the characterization of the efficiency global 
process, considering the assay conditions and reagents all together. It is 
therefore useful for assay optimization once n and k have been fixed and 
to detect changes in rr when k is kept constant or, inversely, changes in k 
while rr is kept unchanged. 

The heterogeneity index rj 

rj measures the internal heterogeneity of the reaction process under 
study. Complex processes include a huge chain of individual and causally 
events inside a multidimensional network of interrelated and interregulated 
biological reactions. Thus, the constant of resistance (k) for the particular 
reaction process under study is a macroscopic indicator of the global 
resistance of that process (k = a1/c1 +a2*2 + ...an/cn/n). If the 
contribution of the individual microscopic constants of resistance (a1/c1, 
a2/c2, ...an/cn) for the individual steps involved in the process were 
homogeneous and no thresholds were present from one step to the next, 
then , no discontinuities in the increase of the Hill coefficient (i.e. in the 
change of k) with R1 should be observed. The existence of a major 
heterogeneity among the k\ values corresponding to the microscopic 
individual steps (i.e. the existence of thresholds for the intermediate 
steps) might lead to a macroscopic discontinuity in the system. 
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Heterogeneity would cause a change in the rate of variation of the Hill 
coefficient and, which would require a jump in the macroscopic value of k 
in order for equation (2) to fit the data. 

The presence of internal heterogeneity in the reaction process can 
5 be detected by the appearance of steps in the rate of change of the Hill 
coefficient, corresponding to the Hill curve that fits the experimental data. 
/7_is defined as an index of heterogeneity and its value corresponds to the 
number of steps in the rate of variation of the Hill coefficient (one step, 
/7=1; two steps, /7 = 2; n steps, /7= n). 
10 o is a key parameter for the dissection and detailed analysis of the 

reaction process. It is useful for the independent optimization and 
development of every one of the steps identified by rj. 

As mentioned, the presence of steps in the rate of change of n 
translates in an abrupt discontinuity in k. Therefore, every step is 
15 determined by a different macroscopic constant of resistance k. Systems 
with /7 = 2, can thus be described by a Hill equation in which k takes two 
different values (*i and k2), according to the R1 interval considered. One 
part of the curve is described by k1 and the other by k2. 

Hill curves describing reaction processes characterized by /7 = 2, are 
20 hybrids generated from two parallel Hill curves differing only in k. The 

transition from one curve to the other may alter the smooth change in the 
slope of the resulting Hill curve. 

The apparent titer r 
In the Hill equation (2), when R1 increases, r increases from 1 to 
25 2,3,4... and P approaches its P max value. On the other direction, on the 
contrary, R1 can only decrease up to a minimal point (R1 min )/ at which r 
and P reach their minimal values. The Hill sigmoidal curve is not 
symmetric, only the right arm is asymptotic (towards P max ). On the left 
arm, the curve has an origin at R1 min ; the empirical curve does not fit the 
30 data for values below R1 
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From a biological point of view, the fact that P does not exist for 
R1 below R1 min , means that there is no 'product' when the concentration 
of 'substrate' is lower than R1 min ; e.g. that the system is not responsive 
to concentrations below R1 min - The minimal concentration of R1 that the 
system can detect and report is R1 min . 

In terms of biological agents, R1 min represents the minimal amount 
that can elicit a response in a given reporter system, and it is represented 
by t. The titer defined this way, is neither an asymptote value nor a value 
approached by extrapolation, but a precise parameter of the Hill equation, 
at the very mathematical origin of the curve. 

t measures the limiting dilution or apparent titer of the biological 
reagent. The value of r is determined by the limit of sensitivity of the 
biological assay system and of method used for the measurement of the 
product P; that is why it is said to be apparent titer. 

t is specific to the batch or stock of the biological reagent tested, r 
represents the apparent concentration of the biological agent and is 
expressed in units per volume, e.g. the maximal dilution of the biological 
agent that leads to the production of P. r is given by the maximal R1 for 
which the Hill coefficient reaches its minimal value (the Hill coefficient 
becomes constant at a value equal or close to 1). The concept of r 
corresponds to that of titer, of general use for viruses, antibodies and 
vectors. Variations in r affect equation (2) by shifting the curve to the 
right or to the left, according to whether the value of r decreases or 
increases, respectively. 

t is a key parameter that measures the 'apparent' concentration of 
a stock of the biological agent, which is necessary for whatever use it will 
be given. 
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Th absolut tit r 0 

9 is a the parameter that measures the absolute concentration 
(titer) of a stock or batch of the biological agent. The value of 9 is not 
determined by nor dependent on the limit of sensitivity of the biological 
assay system or of the method used for the measurement of the product 
P; that is why it is said to be absolute titer. 9 is specific to the batch or 
stock of the biological reagent tested. It represents the real physical 
concentration of the biological agent and is expressed in units per volume, 
e.g. the maximal dilution of the biological agent that leads to the 
production of P. 

9 is given by the following equation 
0/7= r/s (3), 
where s is the sensitivity of the detection method. Therefore, for agents 
detected using the same method, the following expression is valid: 
01/71 /r1 = 92tt2 It2 = 0n/7n/rn = constant (4) 

Using equation (4), the ratio of the absolute titer 9, corresponding 
to two biological agent preparations, can be obtained from their 
respective n and t. Variations in 9 affect the equation (2) by shifting the 
curve to the right or to the left and/or by changing its slope. 

Compensation between rr and #r 

rr and k may appear to compensate to generate two different Hill 
curves (one differing in rr and the other one differing in k) that would 
apparently fit with the same experimental data. As rr and k have opposite 
effects, two Hill curves; in which the increase in rr is compensated by the 
decrease in /c, and vice versa, may seem to represent the same curve, 
which could make it difficult to determine whether two Hill curves are 
different because a change in rr or in k. 

Detailed analysis of the Hill curves indicates that rr and k do not 
compensate very well. Although curves differing in compensatory values 
of either rr or k may vary close each other, they do not fit exactly in any 
of the two regions of highest curvature (before and after the inflection 
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point). This dispersion is caused by the fact that n, but not k, changes 
the slope at the inflection point of the Hill curve. Therefore, e, which is 
the slope of the Hill curve at the inflection point, can be used to easily 
differentiate between two Hill curves that apparently compensate for n 
and k. 

Conclusions 

The application of the Hill analysis to resolve complex biological 
processes is effective for the precise and objective understanding of 
processes like virus or vaccine action, entry, genome replication, 
transgene expression, vector/transgene immunogenicity, cytotoxicity and 
other such parameters. The analysis is independent of the virus vaccine, 
vector and protein type involved and from the output parameter and 
variable measured, such as the internalized vector DNA, transgene mRNA 
level and transgene product activity. 

As in the field of chemical pharmaceuticals, the structure of the 
potential drug (in this case the biological agent) must be optimized to a 
maximal possible intrinsic potency. In analytical development, the goal is 
to search for better performing reporter systems (the lowest possible /c), 
as analytical tool. Two different systems characterized by constants kA 
and /cB, respectively, can be compared (using the same biological agent) 
for their relative resistance or performance. 

Complex systems involving the interaction of biological agents, 
such as viruses, vaccines, gene transfer vectors, antibodies proteins and 
living cells (either in vitro or in vivo) can be analyzed using the Hill 
equation. A complex succession of unitary processes, each of them 
susceptible to be individually analyzed by the Hill equation, as a global 
process, can be also described by the same equation as its constitutive 
steps. 
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EXAMPLE 4 

Materials and Methods 
Cells: 

293 human embryo kidney (HEK) cells, obtained from ATCC, were 
cultured in Dulbecco's modified Eagle's medium containing 4.5 g/l 
glucose (DMEM; GIBGO-BRL) 10 % fetal bovine serum (FBS, Hyclone). 
Hela rep-cap 32 cells, described above, were obtained from Anna Salvetti 
(CHU, Nantes) and cultured in the medium described above. 

Plasmids: 

pNB-Adeno, which encodes the entire E2A and E4 regions and VA 
RNA I and II genes of Adenovirus type 5, was constructed by ligating into 
the polylinker of multiple cloning site of pBSII KS ( + /-) (Stratagene, San 
Diego, USA) the Sall-Hindlll fragment (9842-1 1555 nt) of Adenovirus 
type 5) and the BamHI-Clal fragment (21563- 35950) of pBR325. All 
fragments of adenovirus gene were obtained from the plasmid pBHG-10 
(Microbix, Ontario, Canada). pNB-AAV encodes the genes rep and cap of 
AAV-2 was constructing by ligation of Xbal-Xbal PCR fragment 
containing the genome of AAV-2 from nucleotide 200 to 4480 into Xbal 
site of polylinker MCS of pBSIIKS( + /-). The PCR fragment was obtained 
from pAV1 (ATCC, USA). Plasmid pNB-AAV was derived from plasmid 
pVA1l, which contains the AAV genomic region, rep and cap. pNB-AAV 
does not contain the AAV ITR's present in pAV1. pAAV-CMV(nls)LacZ 
was provided by Dr Anna Salvetti (CHU, Nantes). 

pCMV(nls)LacZ (rAAV vector plasmid) and pNB-Adeno were 
prepared on DH5a E. coli and purified by Nucleobond AX PC500 Kit 
(Macherey-Nagel), according to standard procedures. Plasmid pAAV- 
CMV(nls)LacZ is derived from plasmid psub201 by deleting the rep-cap 
region with SnaB I and replacing it with an expression cassette harboring 
the cytomegalovirus (CMV) immediate early promoter (407 bp), the 
nuclear localized /?-galactosidase gene and the bovine growth hormone 



-79- 



polyA signal (324 bp) (see, Chadeuf eta/. (2000) J. Gene Med. 2:260- 
268. pAAV-CMV(nls)LacZ was provided by Dr Anna Salvetti. 
Virus: 

Wild type adenovirus (AV) type 5 stock, originally provided by Dr 
Philippe Moullier (CHU, Nantes), was produced accordingly to standard 

procedures. 

Construction of Rep mutant libraries 

25 pmol of each mutagenic primer was placed into a 96 PCR well 
plate. 15 a/I of reaction mix (0.25 pmol of pNB-AAV), 25 pmol of the 
selection primer (changing one non-essential unique restriction site to a 
new restriction site), 2 /j\ of 10 X mutagenesis buffer (100 mM Tris- 
acetate pH7.5, 100 mM MgOAc and 500 mM KOAc pH7.5) was added 
into each well. The samples were incubated at 98°C for 5 minutes and 
then immediately incubated for 5 minutes on ice. Finally, the plate was 
placed at room temperature for 30 minutes. 

The primer extension and ligation reactions of the new strands 
were completed by adding to each sample: 7 jj\ of nucleotide mix (2.86 
mM each nucleotide and 1 .43 X mutagenesis buffer) and 3/j\ of a fresh 
1:10 enzyme dilution mix (0.025 U//yl of native T7 DNA polymerase and 
1 U//yl of T4 DNA ligase were diluted in 20 mM Tris HCI pH7.5, 10 mM 
KCI, 10 mM mercaptoethanol, 1 mM DTT, 0.1 mM EDTA and 50% 
glycerol). Samples were incubated at 37°C for 1 hour. The T4 DNA ligase 
was inactivated by incubating the reactions at 72°C for 15 minutes to 
prevent re-ligation of the digested strands during the digestion of the 
parental plasmid (pNB-AAV). 

Each mutagenesis reaction was digested with restriction enzyme to 
eliminate parental plasmids: 30 fj\ solution containing 3 jj\ of 10X enzyme 
digestion buffer and 10 units of restriction enzyme were added to each 
mutagenesis reaction and incubated at 37°C for at least 3 hours. 

90 fj\ of the E. coli XLmutS competent cells (Stratagene, San Diego 
CA; supplemented with 1 .5 //I of jff-mercaptoethanol to a final concen- 
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tration of 25 mM) were aliquoted into prechilled deep-well plates. The 
plates were incubated on ice for 10 minutes and swirling gently every 2 
minutes. 

A fraction of the reactions that had been digested with restriction 
enzyme (1/10 of the total volume) was added to the deep well plates. The 
plates were swirled gently prior to incubation on ice for 30 minutes. A 
heat pulse was performed in a 42°C water bath for 45 seconds, the 
transformation mixture was incubated on ice for 2 minutes and 0.45 ml of 
preheated SOC medium (2% (w/v) tryptone, 0.5% (w/v) yeast extract, 
8.5 mM NaCI, 2.5 mM KCI, 10 mM MgCI 2 and 20 mM glucose at pH 7) 
was added. The plates were incubated at 37°C for 1 hour with shaking. 

To enrich for mutant plasmids, 1 ml of 2X YT broth medium (YT 
medium is 0.5% yeast extract, 0.5% NaCI, 0.8% bacto-tryptone), 
supplemented with 100 jug/ml of ampicillin, was added to each 
transformation mixture and the cultures were grown overnight at 37 °C 
with shaking. Plasmid DNA isolation was performed from each mutant 
culture using standard procedure described in Nucleospin Multi-96 Plus 
Plasmid Kit (Macherey-Nagel). Five hundred /jq of the resulting isolated 
DNA was digested with 10 units of the selection restriction enzyme in a 
total volume of 30/vl containing 3 jj\ of 10X enzyme digestion buffer for 

overnight at 37°C. 

A fraction of the digested reactions (1/10 of the total volume) were 
transformed into 40 /j\ of Epicurian coli XL1-Blue competent cells 
supplemented with 0.68 fj\ of £-mercaptoethanol to a final concentration 
of 25 mM. After heat pulse, 0.45 ml of SOC was added and the 
transformation mixtures were incubated for 1 hour at 37 °C with shaking 
before to be plate on LB-ampicillin agar plates. The agar plates were 
incubated overnight at 37°C and the colonies obtained were picked up 
and grown overnight at 37°C into deep-well plates. 

Four clones per reaction were screened for the presence of the 
mutation using restriction enzyme specific to the new restriction site 
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introduced into the mutated plasmid with the selection primer. The cDNA 
from selected clones was also sequenced to confirm the presence of the 
expected mutation. 

Monitoring rAAV Production 
5 rAAV from each of the above wells, were produced by triple 

transfection on 293 HEK cells. 3 x 10 4 cells were seeded into each well 
of 96 micro-well plate and cultured for 24 hours before transfection. 
Transfection was made on cells at about 70% confluency. 25 kDa PEI 
(poly-ethylene-imine, Sigma-Aldrich) was used for the triple transfection 
10 step. Equimolar amounts of the three plasmids AV helper plasmid (pNB- 



fU Adeno), AAV helper plasmid (pNB-AAV or a mutant clone rep plasmid) 

fi\ and vector plasmid (pAAV-CMV(nls)LacZ) were mixed with 10 mM PEI by 



gently shaking. The mixture was the added to the medium culture on the 
cells. 60 hours after transfection, the culture medium was replaced with 

15 100/vl of lysis buffer (50 mM Hepes, pH 7.4; 150 mM NaCI; 1 m MgCI 2 ; 
1 mM CaCI 2 ; 0.01 % CHAPS). After one cycle of freeze-thawing the 
cellular lysate was filtered through a millipore filter 96 well plate and 
stored at -80°C. 

rAAV infection particles (ip) 

20 Titers of rAAV vector particles were determined on HeLa rep/cap 

32 cells using standard dRA (serial dilution replication assay) test. Cells 
were plated 24 hours before infection at a density of 1 x 10 4 cells in 96- 
well plates. Serial dilutions of the rAAV preparation were made between 1 
and 1 x 10 6 /j\ and used for co-infection of the HeLa rep/cap 32 cells 

25 together with wt-AV type 5 (MOI 25). 48 hours after infection the ip 
were measured by real time PCR or by the quantification of biological 
activity of the transgene. 
Real Time PCR 

Infected HeLa rep/cap 32 cells were lysed with 50 //I of solution 
30 (50 mM Hepes, pH 7.4; 150 mM NaCI). After one cycle of freeze-thawing 
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50 jj\ of Proteinase K (10 mg/ml) and the lysate were incubated one hour 
at 55 °C. The enzyme was inactivated by incubation 10 min at 96°C. 

For real time PCR, 0.2 fj\ of lysate was taken. Final volume of the 
reaction was 10/yl in 384 well plate using an Applied Biosystem Prism 
7900. The primers and fluorescence probe set corresponding to the CMV 
promoter were as follows: CMV 1 primer 5'- 

TGCCAAGTACGCCCCCTAT-3' (SEQ ID No. 733) (0.2 /vM) and CMV 2 
primer 5'-AGGTCATGTACTGGGCATAATGC -3' (SEQ ID No. 734) (0.2 
fjM) ; probe VIC-Tamra 5'-TCAATGACGGTAAATGGCCCGCCT-3' (SEQ 
ID No. 735) (0.1 /jM). dRA plots were obtained by plotting the DNA copy 
number (obtained by real time PCR) vs. the dilution of the rAAV 
preparation. 

/?-Galactosidase activity 

After 48 hours of infection, cells were treated with trypsin, and 
100 jj\ of reaction solution (GalScreen Kit, Tropix) was added and 
incubated for one hour at 26 °C. Luminescence was measured in 
NorthStar (Tropix) HTS station. dRA plots were obtained plotting the 
intensity of £-Galatosidase activity vs. the dilution of the rAAV 
preparation. 

Mathematical Model for results analysis: 

Results were analyzed using the Hill equation-based analysis 
(designated NautScan™; see, Patent n° 9915884, 1999, France; 
published as International PCT application No. WO 01/44809 (PCT n° 
PCT/FR00/03503, Dec, 2000, see EXAMPLE above). Briefly, data were 
processed using a Hill equation-based model that allows extraction of key 
feature indicators of performance for each individual mutant. Mutants 
were ranked based on the values of their individual performance and 
those at the top of the ranking list were selected as Leads. 
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Results 

Generation of diversity. 

To identify candidate amino acid (aa) positions on the rep protein 
involved in rep protein activity an Ala-scan was performed on the rep 
sequence. For this, each amino acid in the rep protein sequence was 
replaced with Alanine. To do this sets of rAAV that encode mutant rep 
proteins in which each differs from wild type by replacement of one 
amino acid with Ala, was generated. Each set of rAAV was individually 
introduced into cells in a well of a microtiter plate, under conditions for 
expression of the rep protein. The amount of virus that could be 
produced from each variant was measured as described below. Briefly, 
activity of Rep was assessed by determining the amount of AAV or rAAV 
produced using infection assays on HeLa Rep-cap 32 cells and by 
measurement of AAV DNA replication using Real Time PCR, or by 
assessing transgene (/?-galactosidase) expression. The relative activity of 
each individual mutant compared to the native protein was assessed and 
"hits" identified. Hit positions are the positions in the mutant proteins 
that resulted in an alteration (selected to be at least about 20%), in this 
instance all resulted in a decrease, in the amount of virus produced 
compared to the activity of the native (wildtype) gene (see Fig. 3A). 

The hits were then used for identification of leads (see, Fig. 3B). 
Assays for Rep activity were performed as described for identification of 
the hit positions. Hit positions on Rep proteins and the effect of specific 
amino acids on the productivity of AAV-2 summarized in the following 
table: 
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| Hit position | 


replacing amino acid (effect) || 


1 4 (ttt) F 1 


(get) A (decrease) | 


1 


1 10 (aag) K I 


(gcg) A (decrease) 




I 20 (ccc) P 


(gec) A (decrease) 


1 


I 22 (att) I I 


(get) A (decrease) 


J 


II 28 (tgg) W I 


(gcg) A (decrease) 




I 32 (gag) E I 


(gcg) A (decrease) 




II 38 (ccg) P 


(gcg) A (decrease) 

* *Lf* ' I 




| 39 (cca) P 


(gca) A (decrease) 


J 


I 54 (ctg) L | 


(get) A (decrease) 


1 


I 59 (ctg) L | 


(gcg) A (decrease) 




I 64 (ctg) L | 


(gcg) A (decrease) 


1 


J 74 (ccg) P I 


(gcg) A (decrease) 


|| 


86 (gag) E 


(gcg) A (decrease) 


I II 


I 88 (tac) Y 


(gec) A (decease) 


1 


I 101 (aaa) K 


(qca) A (decrease) 


1 


I 1 24 (ate) I 


(qcc) A (decrease) 


1 


I 125 (gag) E 


(qcg) A (decrease) 


1 1 


I 1 27 (act) T 


(get) A (decrease) 


1 1 


I 1 32 (ttc) F 


(gec) A (decrease) 


1 1 


I 1 40 (ggc) G 


(gee) A (decrease) 


| J 


I 161 (acc) T 


(gee) A (decrease) 


1 1 


I 1 63 (cct) P 


(get) A (decrease) 


1 1 


I 175 (tat) Y 


(get) A (decrease) 


1 II 


I 193 (ctg) L 


| (gcg) A (decrease) 




I 1 96 (gtg) V 


| (gcg) A (decrease) 




I 1 97 (teg) S 


| (gee) A (decrease) 




I 221 (tea) S 


| (gca) A (decrease) 




I 228 (gtc) V 


| (gcg) A (decrease) 
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Hit position 



replacing amino acid (effect) 



231 (etc) L 



(gec) A (decrease) 



234 (aag) K (gcg) A (decrease) 



237 (acc) T 



(gec) A (decrease) 



250 (tac) Y 



(gec) A (decrease) 



258 (aac) N 



(gec) A (decrease) 



260 (egg) R 



(gcg) A (decrease) 



263 (ate) I 



264 (aag) K 



(gec) A (decrease) 



(gcg) A (decrease) 



334 (ggg) G (gcg) A (decrease) 



335 (cct) V 



(get) A (decrease) 



337 (act) T 



(get) A (decrease) 



341 (acc) T 



(gec) A (decrease) 



342 (aac) N (gec) A (decrease) 



347 (ata) I 



(gca) A (decrease) 



350 (act) T 



(get) A (decrease 



(aat) N (increase) 



354 (tac) Y 



(gec) A (decrease) 



363 (aac) N 



(gec) A (decrease) 



364 (ttt) F 



(get) A (decrease) 



367 (aac) N 



(gee) A (decrease) 



370 (gtc) V 



(gee) A (decrease) 



376 (tgg) W (gcg) A (decrease) 



381 (aag) K (gcg) A (decrease) 



382 (atg) M 



(gcg) A (decrease) 



389 (teg) S 



(gcg) A (decrease) 



407 (tec) S 



(gec) A (decrease) 



41 1 (ata) I 



(gca) A (decrease) 



414 (act) T 



(get) A (decrease) 



420 (tec) S 



(get) A (decrease) 
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Hit position 


replacing amino acid (effect) 


421 (aac) N 


(gec) A (decrease) 




422 (acc) T 


(gec) A (decrease) 




424 (atg) M 


(gcg) A (decrease) 




428 (att) I 


(get) A (decrease) 




429 (gac) D 


(gec) A (decrease) 




438 (cag) Q 


(gcg) A (decrease) 




440 (ccg) P 


(gcg) A (decrease) 




451 (acc) T 


(gec) A (decrease) 




460 (aag) K 


(gcg) A (decrease) 




462 (acc) T 


(gec) A (decrease) 


(ata) I (increase) 


484 (ttc) F 


(gec) A (decrease) 




488 (aag) K 


(gcg) A (decrease) 




495 (ccc) P 


(gec) A (decrease) 




497 (ccc) P 


(gec) A (decrease) 


(cga) R (increase) 


497 (ccc) P 


(gec) A (decrease) 


(etc) L (increase) 


497 (ccc) P 


(gec) A (decrease) 


(tac) Y (increase) 


498 (agt) S 


(get) A (decrease) 




499 (gac) D 


(gec) A (decrease) 




503 (agt) S 


(gcg) A (decrease) 




511 (tea) S 


(gca) A (decrease) 




5 1 2 (gtt) V 


(get) A (decrease) 




516 (teg) S 


(gcg) A (decrease) 




517 (acg) T 


(get) A (decrease) 


(aac) N (increase) 


518 (tea) S 


(gca) A (decrease) 




519 (gac) D 


(gcg) A (decrease) 




542 (ctg) L 


(gcg) A (decrease) 


(teg) S (increase) 


548 (aga) R 


(gca) A (decrease) 


(age) S (increase) 


598 (gga) G 


(gca) A (decrease) 


(gac) D (increase) 
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Hit position 


replacing amino acid (effect) 


598 (gga) G 


(gca) A (decrease) 


(age) S (increase) 




(gcg) A (decrease) 


(ccg) P (increase) 


601 (cca) P 


(gca) A (decrease) 




Hit position 
(within intron) 


replacing sequence (effect) 


630 (tgc) 


gcg (decrease) 


cgc or tea or cct 
(increase) 



The hits in other AAV serotypes (see, also FIG. 4) are as follows 
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HIT POSITION 



/~\/~\ v 


AAV-1 


AAV-3 


AAV-3 B 


AAV-4 


AAV-6 


AAV-5 


A 
*+ 


A 


4 


4 


4 

• 


4 


4 


1 n 


1 0 


10 


10 


10 


10 


10 


9P> 


90 


20 


20 


20 


20 


20 


99 


99 


22 


22 


22 


22 


22 


9Q 


9Q 


29 


29 


29 


29 


29 


32 


32 


32 


32 


32 


32 


32 


38 


38 


38 


38 


38 


38 


38 


39 


39 


39 


39 


39 


39 


39 


54 


54 


54 


54 


54 


54 


54 


59 


59 


59 


59 


59 


59 


59 


64 


64 


64 


64 


64 


64 


64 


74 


74 


74 


74 


74 


74 




86 


86 


86 


86 


86 


86 


85 


88 


88 


88 


88 


88 


88 


87 


101 


101 


101 


101 


101 


101 


100 


124 


124 


124 


124 


124 


124 


123 


125 


125 


125 


125 


125 


125 


124 


127 


127 


127 


127 


127 


127 


126 


132 


132 


132 


132 


132 


132 


131 


140 


140 


140 


140 


140 


140 
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HIT POSITION 


AAV/ O 

AAV-z 


A A \ / 1 

AAV- I 


A A W Q 

AA V-o 


AAV-od 


A A\/ A 


A A\/ 
AMV-D 


AAW C 
AAV-O 


161 


1 61 


161 


I ol 


I bl 


1 Dl 


1 CO 


163 


163 


1 63 


1 63 


-too 

1 63 


1 63 


1 OU 


175 


175 


175 


1 75 


1 75 


1 75 


1 ~7 O 


193 


193 


193 


1 93 


1 93 


■i 00 

1 93 


1 90 


196 


196 


196 


196 


1 96 


1 96 


1 o 

1 93 


197 


197 


197 


197 


1 97 


197 


1 94 


221 


221 


221 


221 


221 


221 


21 7 


228 


228 


228 


O O 

228 


/~i O f""> 

228 


o o o 

228 


o o >i 

224 


231 


231 


231 


231 


231 


O O -1 

231 


O O ~7 

227 


234 


234 


234 


O O >l 

234 


234 


234 


o or\ 
230 


237 


^ 

237 


237 


O O "1 

237 


237 


O 0"7 

237 


o o o 
233 


250 


250 


250 


250 


i it r\ 

250 


O C f\ 

250 


246 


258 


^ ^ 

258 


258 


258 


258 


o c o 

258 


O C A 

254 


-r^v .A 

260 


260 


260 


260 


260 


260 


o c o 

256 


263 


263 


263 


263 


263 


263 


259 


264 


264 


264 


264 


264 


264 


260 


334 


334 


334 


334 


334 


334 


330 


335 


335 


335 


335 


335 


335 


331 


337 


337 


337 


337 


337 


337 


333 


341 


341 


341 


341 


341 


341 


O O ~T 

337 


342 


342 


342 


342 


342 


342 


o o o 

338 


347 


347 


347 


347 


347 


347 


O VI o 

342 


350 


350 


350 


o r" 

350 


350 


350 


346 


354 


354 


354 


354 


354 


O IT VI 

354 


350 


363 


363 


363 


363 


o o o 

363 


oo o 

363 


359 


Jo4 




oD4 




OCA 






367 


367 


367 


367 


367 


367 


363 


370 


370 


370 


370 


370 


370 


366 


376 


376 


376 


376 


376 


376 


372 


381 


381 


381 


381 


381 


381 


377 
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15 



20 



25 
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HIT POSITION 


A A \ / O 

AAV-2 


AAV-1 


AAV-o 


A A \ / OD 
AAV-OD 


A A W A 

AAV-4 


AAV-D 


A AW C 


382 


382 


ft ft ft 

382 


O ft o 

382 


382 


o o o 
OOZ 


o /o 


389 


389 


ft ft r\ 

389 


0 ft ft 
389 


ft O ft 

389 


ft oo 

389 


ft ft rr 
OOD 


407 


407 


407 


407 


A ft ~~l 

407 


A ft ~7 

407 


40o 


41 1 


41 1 


41 1 


41 1 


41 4 4 

41 1 


i 4 4 

41 1 


407 


414 


414 


41 A A 

414 


414 


AAA 

414 


414 


410 


420 


420 


420 


420 


A ft ft 

420 


/■ o ft 
420 


41 6 


421 


421 


421 


A A 

421 


A ft <■ 

421 


421 


41 7 


422 


422 


422 


422 


ji ft ft 
422 


41 ft ft 

422 


41 8 


424 


424 


A f~\ A 

424 


424 


^ ft >i 
424 


>1 ft yl 

424 


A ft ft 

420 


428 


428 


428 


428 


>i ft ft 
428 


/•Oft 

428 


A ft /I 

424 


429 


429 


429 


A ft ft 

429 


A ft ft 

429 


>l ft ft 

429 


4zo 


438 


438 


438 


A ft ft 

438 


>i ft ft 
438 


A ft ft 

438 


4o4 


440 


440 


440 


440 


440 


440 


A ft 

436 


451 


451 


451 


451 


451 


451 


44—7 
447 


460 


460 


460 


460 


j* ft^\ 
460 


A ft ^\ 

460 


4 r 0 

456 


462 


462 


^ *»V 4**\ 

462 


462 


462 


462 


41 f~ ft 

458 | 


484 


484 


484 


484 


484 


Jt ft 41 

484 


41 ftft 

480 


488 


488 


488 


488 


488 


488 


41 ft 41 

484 


495 


495 


495 


495 


495 


ji ^% f~ 

495 


A f\ A 

491 


497 


497 


497 


497 


ji r\ 

497 


41 ft — 1 

497 


41 ft ft 

493 


498 


498 


498 


498 


498 


A ft ft 

498 


A ft A 

494 


499 


M 4*^ A 

499 


499 


499 


499 


A ft ft 

499 


a ft tr 

495 


503 


503 


503 


503 


503 


IT ft ft 

503 


499 


51 1 


51 1 


51 1 


|— A A 

51 1 


51 1 


P 4 4 

511 


C ft ft 

529 


512 


512 


512 


^ A ft 

51 2 


^ ft 
51 2 


C 1 ft 

512 


er ft ft 

530 


516 


DID 


bio 


DID 


D 1 D 


0 1 O 


DO 4 !- 


517 


517 


517 


517 


517 


517 


535 


518 


518 


518 


518 


518 


518 


536 


519 


519 


519 


519 


519 


519 


537 


542 


543 


542 


542 


542 


543 


561 
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HIT POSITION 


AAV-2 


A A \ # 1 

AAV-1 


AAV/ O 

AAV-3 


A A \ / o 

AAV-3B 


A A \ / A 

AAV-4 


A A \ / /"> 

AAV-6 


AAV/ r - 

AAV-5 


548 


549 


548 


548 


548 


549 


567 


598 


599 


600 


600 


599 


599 




600 


602 


603 


603 


602 


602 


589 


601 


603 


604 


604 


603 


603 


590 



f \ s 

5 "K 1 

a.- 

-i 



Sets of nucleic acids encoding the rep protein were generated. The 
rep proteins encoded by these sets of nucleic acid molecules were those 
in which each amino acid position identied as a "hit" in the ala-scan step, 
were each sequentially replaced by all remaining 1 8 amino acids using 
10 site directed mutagenesis. Each mutant was designed, generated, 
processed and analyzed physically separated from the others in 
addressable arrays. No mixtures, pools, nor combinatorial processing 
were used. 

As in the first round (alanine scan), a library of mutant rAAV was 
15 generated in which each individual mutant was independently and 
individually generated in a independent reaction and such that each 
mutant contains only a single amino acid change and this for each amino 
acid residue. Again, each resulting mutant rep protein was then 
expressed and the amount of virus produced in cells assessed and 
20 compared to the native protein. 
Lead identification 

» 

Since rep proteins that result in increased virus production are of 
interest, those mutants that lead to an increase in the amount of virus 
produced (2 to 10 times the native activity), were selected as "leads." 
25 Ten such mutants were identified. 

Based on the results obtained from the assays described above (i.e. 
titer of virus produced by each rep variant), each individual rep variant 
was assigned a specific activity. Those variant proteins displaying the 
highest titers were selected as leads (see Table above). Leads include: 



-91 



37851-911 



amino acid replacement of T by N at Hit position 350; T by I at Hit 
position 462; P by R at Hit position 497; P by L at Hit position 497; P by 
Y at Hit position 497; T by N at Hit position 517; G by S at Hit position 
598; G by D at Hit position 598; V by P at Hit position 600. 
5 Also provided are combinations of the above mutant Rep 78, 68, 

52. 40 proteins, nucleic acids encoding the proteins, and recombinant 
AAV (any serotype) contains the mutation at the indicated position or 
corresponding position for serotypes other than AAV-2, including any set 
forth in the following table and corresponding SEQ ID Nos. Each amino 
O 10 acid sequence is set forth in a separate sequence ID listing; for each 

mutation or combination thereof there is a single SEQ ID setting forth the 
unspliced nucleic acid sequence for Rep78/68, which for all mutations 
from amino acid 228 on, includes the corresponding Rep 52 and Rep 40 
encoding sequence as well. 
1 5 Amino acid sequences of exemplary mutant Rep proteins 
Seq no. gene position(s) codon(s) 



rS ! 

5 b : 

ill 



■as. 
- e t 

t 

K23 



a s ; 



i s 

£5 S3 



5 ! 



seq.1 rep78 4 GCT 

seq. 2 rep68 4 GCT 

seq. 3 rep78 10 GCG 

20 seq. 4 rep68 10 GCG 

seq. 5 rep78 20 GCC 

seq. 6 rep68 20 GCC 

seq. 7 rep78 22 GCT 

seq. 8 rep68 22 GCT 

25 seq. 9 rep78 29 GCG 

seq. 10 rep68 29 GCG 

seq. 11 rep78 38 GCG 

seq. 12 rep68 38 GCG 

seq. 13 rep78 39 GCA 

30 seq. 14 rep68 39 GCA 

seq. 15 rep78 53 GCT 

seq. 16 rep68 53 GCT 

seq. 17 rep78 59 GCG 

seq. 18 rep68 59 GCG 

35 seq. 19 rep78 64 GCT 

seq. 20 rep68 64 GCT 

seq. 21 rep78 74 GCG 

seq. 22 rep68 74 GCG 

seq. 23 rep78 86 GCG 

40 seq. 24 rep68 86 GCG 

seq. 25 rep78 88 GCC 

seq. 26 rep68 88 GCC 

seq. 27 rep78 101 GCA 
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seq.28 rep68 

seq.29 rep78 

seq.30 rep68 

seq.31 rep78 

5 seq.32 rep68 

seq.33 rep78 

seq.34 rep68 

seq.35 rep78 

seq.36 rep68 

10 seq.37 rep78 

seq.38 rep68 

seq.39 rep78 

seq.40 rep68 

seq.41 rep78 

15 seq.42 rep68 

seq.43 rep78 

seq.44 rep68 

seq.45 rep78 
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4 O ~T JI O O 

1 97 420 


ooo ooo 

GCG GCC 


seq.452 


jo\ 

rep40 


<1 O T JI o o 

1 97 420 


ooo ooo 

GCG GCC 


seq.453 


rep78 


f— m OOO /IOC 

54 338 495 


ooo ooo ooo 

GCC GCC GCC 


seq.454 


rep52 


b4 ooo 4yo 


rrr rrr rrr 

OLL VJlrV^ JjV^L* 


seq.455 


rep68 


54 338 495 


GCC GCC GCC 


seq.456 


rep40 


54 338 495 


GCC GCC GCC 


seq.457 


rep78 


197 427 


GCG GCG 


seq.458 


rep 5 2 


197 427 


GCG GCG 


seq.459 


rep68 


197 427 


GCG GCG 
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:. : 
fas' 

f\ \ 

3 5 i 

5 s ? I 
3 fe' 
pi 

*■ 

. («! 
: 



1 U 



7! 



:=Sf 



seq.460 
seq.461 
seq.462 
seq.463 
5 seq.464 
seq.465 
seq.466 
seq.467 
seq.468 

10 seq.469 
seq.470 
seq.471 
seq.472 
seq.473 

15 seq.474 
seq.475 
seq.476 
seq.477 
seq.478 

20 seq.479 
seq.480 
seq.481 
seq.482 
seq.483 

25 seq.484 
seq.485 
seq.486 
seq.487 
seq.488 

30 seq.489 
seq.490 
seq.491 
seq.492 
seq.493 

35 seq.494 
seq.495 
seq.496 
seq.497 
seq.498 

40 seq.499 
seq.500 
seq.501 
seq.502 
seq.503 

45 seq.504 
seq.505 
seq.506 
seq.507 
seq.508 

50 seq.509 
seq.510 
seq.51 1 
seq.51 2 
seq.513 



rep40 

rep78 

rep52 

rep68 

rep40 

rep78 

rep52 

rep68 

rep40 

rep78 

rep68 

rep78 

rep52 

rep68 

rep40 

rep78 

rep52 

rep68 

rep40 

rep78 

rep52 

rep68 

rep40 

rep78 

rep68 

rep78 

rep52 

rep68 

rep40 

rep78 

rep52 

rep68 

rep40 

rep78 

rep52 

rep68 

rep40 

rep78 

rep68 

rep78 

rep52 

rep68 

rep40 

rep78 

rep52 

rep68 

rep40 

rep78 

rep52 

rep68 

rep40 

rep78 

rep52 

rep68 



197 427 

54 228 370 387 

54 228 370 387 

54 228 370 387 

54 228 370 387 

221 289 

221 289 

221 289 

221 289 

54 163 

54 163 

341 407.420 

341 407 420 

341 407 420 

341 407 420 

54 228 

54 228 

54 228 

54 228 

96 125 511 

96 125 511 

96 125 511 

96 125 511 

54 163 

54 163 

1 97 420 

197 420 

197 420 

197 420 

334 428 499 

334 428 499 

334 428 499 

334 428 499 

197 414 

197 414 

197 414 

197 414 

30 54 127 

30 54 127 

29 260 

29 260 

29 260 

29 260 

4 484 

4 484 

4 484 

4 484 

258 124 132 
258 124 132 
258 124 132 
258 124 132 
231 497 
231 497 
231 497 



GCG GCG 

GCC GCC GCC GCG 

GCC GCC GCC GCG 

GCC GCC GCC GCG 

GCC GCC GCC GCG 

GCA GCC 

GCA GCC 

GCA GCC 

GCA GCC 

GCC GCT 

GCC GCT 

GCC GCC GCC 

GCC GCC GCC 

GCC GCC GCC 

GCC GCC GCC 

GCC GCC 

GCC GCC 

GCC GCC 

GCC GCC 

GCA GCG GCA 

GCA GCG GCA 

GCA GCG GCA 

GCA GCG GCA 

GCC GCT 

GCC GCT 

GCG GCC 

GCG GCC 

GCG GCC 

GCG GCC 

GCG GCT GCC 

GCG GCT GCC 

GCG GCT GCC 

GCG GCT GCC 

GCG GCT 

GCG GCT 

GCG GCT 

GCG GCT 

GCG GCC GCT 

GCG GCC GCT 

GCG GCG 

GCG GCG 

GCG GCG 

GCG GCG 

GCT GCC 

GCT GCC 

GCT GCC 

GCT GCC 

GCC GCC GCC 

GCC GCC GCC 

GCC GCC GCC 

GCC GCC GCC 

GCC GCC 

GCC GCC 

GCC GCC 
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seq.514 rep40 

seq.515 rep78 

seq.516 rep52 

seq.517 rep68 

5 seq.518 rep40 

seq.519 rep78 

seq.520 rep52 

seq.521 rep68 

seq.522 rep40 

10 seq.523 rep78 

seq.524 rep52 

seq.525 rep68 

seq.526 rep40 

seq.527 rep78 

15 seq.528 rep68 

M seq.529 rep78 

P seq.530 rep52 

P seq.531 rep68 

PJ seq.532 rep40 

fU 20 seq.533 rep78 

f|j seq.534 rep52 

Jg. seq.535 rep68 

yg seq.536 rep40 

~ seq.537 rep78 

y, 25 seq.538 rep52 

hi seq.539 rep68 

y seq.540 rep40 

[J: seq.541 rep78 

J! seq.542 rep52 

H 30 seq.543 rep68 

r * seq.544 rep40 

seq.545 rep78 

seq.546 rep52 

seq.547 rep68 

35 seq.548 rep40 

seq.549 rep78 

seq.550 rep52 

seq.551 rep68 

seq.552 rep40 

40 seq.553 rep78 

seq.554 rep52 

seq.555 rep68 

seq.556 rep40 

seq.557 rep78 

45 seq.558 rep52 

seq.559 rep68 

seq.560 rep40 

seq.561 rep78 

seq.562 rep68 

50 DNA Sequences 

Sequence aa position 

seq.563 4 

seq.564 10 



231 497 


GCC GCC 


221 258 


GCA GCC 


221 258 


GCA GCC 


221 258 


GCA GCC 


221 258 

■ ■ X>/ X^ 


GCA GCC 


234 264 326 

^— x^ ■ ^— X# ■ X^ ^— X^ 


GCG GCG GCC 


234 264 326 

X^ I Xj^ ¥ Xj^ X^ 


GCG GCG GCC 


234 264 326 

I * Xj^ ■ X^ ^— X^ 


GCG GCG GCC 


234 264 326 

X^ 1 W A— xx ■ x^ 


GCG GCG GCC 


153 398 

1 **** V/ X^ 


AGC GCG 


1 53 398 

1 xV x< xv w 


AGC GCG 

w * xj x^ x^ 


153 398 

■ x** v./ x>/ x/ xv 


AGC GCG 


1 53 398 


AGC GCG 


53 216 

w XV a_ 1 w 


GCG GCC 


53 21 6 


GCG GCC 


22 382 


GCT GCG 

XJ 1 x* x^ x^ 


22 382 


GCT GCG 


22 382 


GCT GCG 

Xm* X^ 1 U w 


22 382 

A— Ak. WW A- 


GCT GCG 


231 41 1 

a_ xv i r i ■ 


GCC GCA 


231 41 1 

£— XV I ■ 1 I 


GCC GCA 


231 41 1 

^— xv 1 ■ 1 V 


GCC GCA 


231 41 1 


GCC GCA 


59 305 

xV v/ x^ xv w 


GCG GCC 


59 305 


GCG GCC 


59 305 

w xv XV xv v/ 


GCG GCC 

XJ X^ > x^ 


59 305 

x/ xv xV xv w 


GCG GCC 

> ■ * Xj^ > »pJ X»^ X^ Xj^ 


53 231 

xV XV ^— XV 1 


GCG GCC 

XJ Xj^ X «* > ^ Xj^ Xj^ 


53 231 

XV XV £— XV 1 


GCG GCC 

X r X/ x»j x^ x^ 


53 231 


GCG GCC 


53 231 


GCG GCC 

XJ X^ XmJ W X* XV 


258 498 


GCC GCT 

X*J VV XjV XJ X»# 1 


258 498 


GCC GCT 

XJ W X^ XJ X»r I 


258 498 


GCC GCT 

xJ xv xv x*^ Xrf I 


258 498 


GCC GCT 

XV XV XV XV w ff 


88 231 


GCC GCC 

XV XV XV XV X^ Xj^ 


88 231 


GCC GCC 

XJ xv xv W XV XV 


88 231 


GCC GCC 

XV XV ^V Xrf X^ 


88 231 


GCC GCC 

xv x^ X*-f x»^ 


101 363 


GCA GCC 

xv x^ ■ i x«^ 


101 363 

i xv i xv xv ^v 


GCA GCC 

X^^ Xj^ # % x«# x^ x^ 


101 363 

1 XV 1 X./ XV XV 


GCA GCC 


101 ^fi^ 

1 V_/ 1 ouo 


fiCA GCC 


354 132 


GCC GCC 


354 132 


GCC GCC 


354 132 


GCC GCC 


354 132 


GCC GCC 


10 132 


GCG GCC 


10 132 


GCG GCC 



codon 

GCT 

GCG 
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10 



15 



20 



Jssa. 

I s 



i Is' 
n ! 



JL. 



: a. 



3 a. 



u 

w 30 

Li. 



25 



35 



45 



50 



seq.565 


20 


GCC 


seq.566 


22 


GCT 


seq.567 


29 


GCG 


seq.568 


38 


GCG 


seq.569 


39 


GCA 


seq.570 


53 


GCT 


seq.571 


59 


f\ f\ 

GCG 


seq.572 


64 


GCT 


seq.573 


74 


GCG 


seq.574 


86 


GCG 


seq.575 


88 


j^-V *~*S 

GCC 


seq.576 


101 


GCA 


seq.577 


124 


GCC 


seq.578 


125 


GCG 


seq.579 


127 


GCT 


seq.580 


132 


GCC 


seq.581 


140 


GCC 


seq.582 


161 


GCC 


seq.583 


163 


GCT 


seq.584 


175 


GCT 


seq.585 


193 


GCG 


seq.586 


196 


GCC 


seq.587 


197 


y^^v y"x 

GCC 


seq.588 


221 


GCA 


seq.589 


228 (Rep78/68) 


GCG 




228 (Rep52) 


GCG 




228 (Rep 40) 


GCG 


seq.590 


231 (Rep78/68) 


GCC 




231 (Rep 52) 


GCC 




231 (Rep 40) 


GCC 


seq.591 


234 (Rep78/68) 


GCG 




234 (Rep 52) 


GCG 




234 (Rep 40) 


GCG 


seq.592 


237 (Rep78/68) 


GCC 




237 (Rep 52) 


GCC 




237 (Rep 40) 


GCC 


seq.593 


250 (Rep78/68) 


GCC 




250 


GCC 




250 


GCC 


seq.594 


258 (Rep78/68) 


GCC 




258 


a y"* 

GCC 




258 


GCC 


seq.595 


260 (Rep78/68) 


GCG 




260 


GCG 




260 


GCG 


seq.596 


263 (Rep78/68) 


GCC 




263 


GCC 




263 


GCC 


seq.597 


264 (Rep78/68) 


GCG 




264 


GCG 




264 


GCG 


seq.598 


334 (Rep78/68) 


GCG 




334 


GCG 




334 


GCG 



-10: 



seq.599 335 (Rep78/68) GCT 

335 GCT 

335 GCT 

seq.600 337 (Rep78/68) GCT 

5 337 GCT 

337 GCT 

seq.601 341 (Rep78/68) GCC 

341 GCC 

341 GCC 
10 seq.602 342 (Rep78/68) GCC 

342 GCC 
342 GCC 

seq.603 347 (Rep78/68) GCA 

347 GCA 

1 5 347 GCA 

seq.604 350 (Rep78/68) AAT 

350 AAT 

350 AAT 

seq.605 350 (Rep78/68) GCT 

20 350 GCT 

350 GCT 

seq.606 354 (Rep78/68) GCC 

354 GCC 

354 GCC 

25 seq.607 363 (Rep78/68) GCC 

363 GCC 

363 GCC 
seq.608 364 (Rep78/68) GCT 

364 GCT 
30 364 GCT 

seq.609 367 (Rep78/68) GCC 

367 GCC 

367 GCC 

seq.610 370 (Rep78/68) GCC 

35 370 GCC 

370 GCC 

seq.611 376 (Rep78/68) GCG 

376 GCG 

376 GCG 

40 seq.612 381 (Rep78/68) GCG 

381 GCG 

381 GCG 
seq.613 382 (Rep78/68) GCG 

382 GCG 
45 382 GCG 

seq.614 389 (Rep78/68) GCG 

389 GCG 

389 GCG 

seq.615 407 (Rep78/68) GCC 

50 407 GCC 

407 GCC 

seq.616 411 (Rep78/68) GCA 

411 GCA 

41 1 GCA 
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5 
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3 Li 



3 



20 



30 



35 



45 



50 



seq.617 


414 (Rep78/68) 


GCT 




414 


GCT 




414 


GCT 


seq.618 


420 (Rep78/68) 


GCT 




420 


GCT 




420 


GCT 


seq.61 9 


421 (Rep78/68) 


GCC 




421 


GCC 




421 


GCC 


seq.620 


422 (Rep78/68) 


GCC 




422 


GCC 




422 


GCC 


seq.621 


424 (Rep78/68) 


GCG 




424 


GCG 




424 


GCG 


seq.622 


428 (Rep78/68) 


GCT 




428 


GCT 




428 


GCT 


seq.623 


429 (Rep78/68) 


GCC 




429 


GCC 




429 


GCC 


seq.624 


438 (Rep78/68) 


GCG 


438 


GCG 




438 


GCG 


seq.625 


440 (Rep78/68) 


GCG 




440 


GCG 




440 


GCG 


seq.626 


451 (Rep78/68) 


GCC 




451 


GCC 




451 


GCC 


seq.627 


460 (Rep78/68) 


GCG 




460 


GCG 




460 


GCG 


seq.628 


462 (Rep78/68) 


GCC 




462 


GCC 




462 


GCC 


seq.629 


462 (Rep78/68) 


ATA 




462 


ATA 




462 


ATA 


seq.630 


484 (Rep78/68) 


GCC 




484 


GCC 




484 


GCC 


seq.631 


488 (Rep78/68) 


GCG 


488 


y*™* y"\ 

GCG 




488 


GCG 


seq.632 


495 <Rep78/68) 


GCC 




495 


GCC 




495 


GCC 


seq.633 


497 (Rep78/68) 


GCC 




497 


GCC 




497 


GCC 


seq.634 


497 (Rep78/68) 


CGA 




497 


CGA 




497 


CGA 



-1 




seq.635 497 (Rep78/68) CTC 

497 CTC 

497 CTC 

seq.636 497 (Rep78/68) TAC 

5 497 TAC 

497 TAC 
seq.637 498 (Rep78/68) GCT 

498 GCT 

498 GCT 
10 seq.638 499 (Rep78/68) GCC 

499 GCC 
499 GCC 

seq.639 503 (Rep78/68) GCG 

503 GCG 

15 503 GCG 

seq.640 510 (Rep78/68) GCA 

D 510 GCA 

D 510 GCA 

HI seq.641 511 (Rep78/68) GCA 

n\ 20 511 GCA 

ni 511 GCA 

j? seq.642 512 (Rep78/68) GCT 

>! 512 GCT 

•Has' 

512 GCT 

L, 25 seq.643 516 (Rep78/68) GCG 

516 GCG 
:'i 516 GCG 

seq.644 517 (Rep78/68) GCT 

517 GCT 
! S J 30 517 GCT 
|s ~ : seq.645 51 7 (Rep78/68) AAC 

517 AAC 

517 AAC 
seq.646 518 (Rep78/68) GCA 

35 518 GCA 

518 GCA 
seq.647 519 (Rep78/68) GCG 

519 GCG 
519 GCG 

40 seq.648 598 (Rep78/68) GCA 

seq.649 600 (Rep78/68) GCG 

seq.650 601 (Rep78/68) GCA 

seq.651 335 420 495 GCT GCC GCC 

335 420 495 GCT GCC GCC 

45 335 420 495 GCT GCC GCC 

seq.652 39 140 GCA GCC 

seq.653 279 428 451 GCC GCT GCC 

279 428 451 GCC GCT GCC 

279 428 451 GCC GCT GCC 

50 seq.654 125 237 600 GCG GCC GCG 

1 25 237 600 GCG GCC GCG 

1 25 237 600 GCG GCC GCG 

seq.655 163 259 GCT GCG 

163 259 GCT GCG 
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163 259 GCTGCG 

seq.656 17 127 189 GCG GCT GCG 

seq.657 350 428 GCT GCT 

350 428 GCT GCT 

5 350 428 GCT GCT 

seq.658 54 338 495 GCC GCC GCC 

54 338 495 GCC GCC GCC 

54 338 495 GCC GCC GCC 

seq.659 350 420 GCT GCC 

10 350 420 GCT GCC 

350 420 GCT GCC 

seq.660 189 197 518 GCG GCG GCA 

189 197 518 GCG GCG GCA 

189 197 518 GCG GCG GCA 

15 seq.661 468 516 GCC GCG 

468 516 GCC GCG 

468 516 GCC GCG 
seq.662 127 221 350 54 140 GCT GCA GCT GCC GCC 

127 221 350 54 140 GCT GCA GCT GCC GCC 
20 1 27 221 350 54 140 GCT GCA GCT GCC GCC 

seq.663 221 285 GCA GCG 

221 285 GCA GCG 

221 285 GCA GCG 

seq.664 23 495 GCT GCC 

25 23 495 GCT GCC 

23 495 GCT GCC 

P seq.665 20 54 420 495 GCC GCC GCC GCC 

20 54 420 495 GCC GCC GCC GCC 

20 54 420 495 GCC GCC GCC GCC 

30 seq.666 412 612 GCC GCG 

412 612 GCC GCG 

412 612 GCC GCG 

seq.667 1 97 41 2 GCG GCC 

197 412 GCG GCC 

35 197 412 GCG GCC 

seq.668 412 495 511 GCC GCC GCA 

41 2 495 51 1 GCC GCC GCA 

41 2 495 51 1 GCC GCC GCA 

seq.669 98 422 GCC GCC 

40 98 422 GCC GCC 

98 422 GCC GCC 

seq.670 1 7 1 27 1 89 GCG GCT GCG 

seq.671 20 54 495 GCC GCC GCC 

20 54 495 GCC GCC GCC 

45 20 54 495 GCC GCC GCC 

seq.672 54 163 GCC GCT 

seq.673 259 54 GCG GCC 

259 54 GCG GCC 

259 54 GCG GCC 

50 seq.674 335 399 GCT GCG 

335 399 GCT GCG 

335 399 GCT GCG 

seq.675 221 432 GCA GCA 

221 432 GCA GCA 
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221 432 


GCA GCA 


seq.676 


259 516 


GCG GCG 




259 516 


GCG GCG 




259 516 


GCG GCG 


seq.677 


495 516 


GCC GCG 




495 516 


GCC GCG 




495 516 


GCC GCG 


seq.678 


414 14 


GCT GCC 




414 14 


GCT GCC 




414 14 


GCT GCC 


seq.679 


74 402 495 


GCG GCC GCC 


74 402 495 


GCG GCC GCC 




74 402 495 


GCG GCC GCC 


seq.680 


228 462 497 


GCC GCC GCC 




228 462 497 


GCC GCC GCC 




228 462 497 


GCC GCC GCC 


seq.681 


290 338 


GCG GCC 




290 338 


GCG GCC 




290 338 


GCG GCC 


seq.682 


140 511 


GCC GCA 




140 511 


GCC GCA 




140 511 


GCC GCA 


seq.683 


86 378 


GCG GCG 


86 378 


GCG GCG 




86 378 


GCG GCG 


seq.684 


54 86 


GCC GCG 




54 86 


GCC GCG 




54 86 


GCC GCG 


seq.685 


214 495 140 


GCG GCC GCC 


214 495 140 


GCG GCC GCC 




214 495 140 


GCG GCC GCC 


seq.686 


495 511 


GCC GCA 




495 511 


GCC GCA 




495 511 


GCC GCA 


seq.687 


495 54 


GCC GCC 




495 54 


GCC GCC 




495 54 


GCC GCC 


seq.688 


197 495 


GCG GCC 




197 495 


GCG GCC 




197 495 


GCG GCC 


seq.689 


261 20 


GCC GCC 




261 20 


GCC GCC 




261 20 


GCC GCC 


seq.690 


54 20 


GCC GCC 


seq.691 


197 420 


GCG GCC 




197 420 


GCG GCC 




197 420 


GCG GCC 


seq.692 


54 338 495 


GCC GCC GCC 




54 338 495 


GCC GCC GCC 




54 338 495 


GCC GCC GCC 


seq.693 


197 427 


GCG GCG 




197 427 


GCG GCG 




197 427 


GCG GCG 


seq.694 


54 228 370 387 


GCC GCC GCC GCG 
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54 228 370 387 


GCC GCC GCC GCG 




54 228 370 387 


GCC GCC GCC GCG 


seq.695 


221 289 


GCA GCC 


221 289 


GCA GCC 




221 289 


GCA GCC 


seq.696 


54 163 


GCC GCT 


54 163 


GCC GCT 


seq.697 


341 407 420 


GCC GCC GCC 


341 407 420 


GCC GCC GCC 




341 407 420 


GCC GCC GCC 


seq.698 


54 228 


GCC GCC 


54 228 


GCC GCC 




54 228 


GCC GCC 


seq.699 


96 125 511 


GCA GCG GCA 


96 125 511 


GCA GCG GCA 




96 125 51 1 


GCA GCG GCA 


seq.700 


1 97 420 


GCG GCC 


197 420 


GCG GCC 




197 420 


GCG GCC 


seq.701 


334 428 499 


GCG GCT GCC 


334 428 499 


GCG GCT GCC 




334 428 499 


GCG GCT GCC 


seq.702 


197 414 


GCG GCT 


197 414 


GCG GCT 




197 414 


GCG GCT 


seq.703 


30 54 1 27 


GCG GCC GCT 


seq.704 


29 260 


GCG GCG 


29 260 


GCG GCG 




29 260 


GCG GCG 


seq.706 


4 484 


GCT GCC 


4 484 


GCT GCC 




4 484 


GCT GCC 


seq.707 


258 124 132 


GCC GCC GCC 


258 124 132 


GCC GCC GCC 




258 124 132 


GCC GCC GCC 


seq.708 


231 497 


GCC GCC 


231 497 


GCC GCC 




231 497 


GCC GCC 


seq.709 


221 258 


GCA GCC 


221 258 


GCA GCC 




221 258 


GCA GCC 


seq.71 0 


234 264 326 


GCG GCG GCC 


234 264 326 


GCG GCG GCC 




234 264 326 


GCG GCG GCC 


seq.71 1 


153 398 


AGC GCG 


153 398 


AGC GCG 




153 398 


AGC GCG 


seq.71 2 


53 216 


GCG GCC 


seq.71 3 


22 382 


GCT GCG 


22 382 


GCT GCG 




22 382 


GCT GCG 


seq.714 


231 411 


GCC GCA 


231 411 


GCC GCA 




231 41 1 


GCC GCA 
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seq.715 


59 305 


GCG GCC 




59 305 


GCG GCC 




59 305 


GCG GCC 


seq.71 6 


53 231 


GCG GCC 




53 231 


GCG GCC 




53 231 


GCG GCC 


seq.71 7 


258 498 


GCC GCT 




258 498 


GCC GCT 




258 498 


GCC GCT 


seq.718 


88 231 


GCC GCC 




88 231 


GCC GCC 




88 231 


GCC GCC 


seq.71 9 


101 363 


GCA GCC 




101 363 


GCA GCC 




101 363 


GCA GCC 


seq.720 


354 132 


GCC GCC 




354 132 


GCC GCC 




354 132 


GCC GCC 


seq.726 


598 


GAC 


seq.727 


598 


AGC 


seq.728 


600 


CCG 



The above nucleic acid molecules are provided in plasmids, which 
are introduced into cells to produce the encoded proteins. The analysis 
revealed the amino acid positions that affect Rep proteins activities. 
Changes of amino acids at any of the hit positions result in altered protein 
activity. Hit positions are numbered and referenced starting from amino 
acid 1 (nucleotide 321 in AAV-2 genome), also codon 1 of the protein 
Rep78 coding sequence under control of p5 promoter of AAV-2: 4, 20, 
22, 29, 32, 38, 39, 54, 59, 124, 125, 127, 132, 140, 161, 163, 193, 
196, 197, 221, 228, 231, 234, 258, 260, 263, 264, 334, 335, 337, 
342, 347, 350, 354, 363, 364, 367, 370, 376, 381, 389, 407, 411, 
414, 420, 421, 422, 424, 428, 438, 440, 451, 460, 462, 484, 488, 
495, 497, 498, 499, 503, 511, 512, 516, 517, 518, 542, 548, 598, 
600 and 601 . The encoded Rep78, Rep68, Rep 52 and Rep 40 proteins 
and rAAV encoding the mutant proteins are provided. The corresponding 
nucleic acid molecules, Rep proteins, rAAV and cells containing the 
nucleic acid molecules or rAAV in which the native proteins are from 
other AAV serotypes, including, but are not limited to, AAV-1, AAV-3, 
AAV-3B, AAV-4, AAV-5 and AAV-6. 
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Other hit positions identified include: 10, 64, 74, 86, 88, 101, 
175, 237, 250, 334, 429 and 519. 

Also provided are nucleic acid molecules, the rAAV, and the 
encoded proteins in which the native amino acid at each hit position is 
replaced with another amino acid, or is deleted, or contains additional 
amino acids at or adjacent to or near the hit positions. In particular the 
following nucleic acid molecules and rAAV that encode proteins 
containing the following amino acid replacements or combinations 
thereof: T by N at Hit position 350; T by I at Hit position 462; P by R at 
Hit position 497; P by L at Hit position 497; P by Y at Hit position 497; T 
by N at Hit position 517; L by S at hit position 542; R by S at hit position 
548; G by D at Hit position 598; G by S at Hit position 598; V by P at Hit 
position 600; in order to increase Rep proteins activities in terms on AAV 
or rAAV productivity. The corresponding nucleic acid molecules, 
recombinant Rep proteins from the other serotypes and the resulting 
rAAV are also provided (see Figs. 5 and the above Table for the 
corresponding position in AAV-1, AAV-3, AAV-3B, AAV-4, AAV-5 and 
AAV-6). 

Mutant adeno-associate virus (AAV) Rep proteins and viruses 
encoding such proteins that include mutations at one or more of residues 
64, 74, 88, 175, 237, 250 and 429, where residue 1 corresponds to 
residue 1 of the Rep78 protein encoding by nucleotides 321-323 of the 
AAV-2 genome, and where the amino acids are replaced as follows: L by 
A at position 64; P by A at position 74; Y by A at position 88; Y by A at 
position 175; T by A at position 237; T by A at position 250; D by A at 
position 429 are provided. Nucleic acid molecules encoding these 
viruses and the mutant proteins are also provided. 

Also provided are nucleic acid molecules produced from any of the 
above-noted nucleic acid molecules by any directed evolution method, 
including, but are not limited to, re-synthesis, mutagenesis, recombination 
and gene shuffling and any way by combining any combination of the 
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molecules, i.e., one, two by one, two by two, n by n, where n is the 

number of molecules to be combined ( i.e., combining all together). The 
resulting recombinant AAV and encoded proteins are also provided. 

Also provided are nucleic acid molecule in which additional amino 
acids surrounding each hit, such as one, two, three . . . ten or more, 
amino acids are systematically replaced, such that the resulting Rep 
protein(s) has increased or decreased activity. Increased activity as 
assessed by increased recombinant virus production in suitable cells is of 
particular interest for production of recombinant viruses for use, for 
example, in gene therapy. 

Also provided are combinations of the above noted mutants in 
which several of the noted amino acids are changed and optionally 
additional amino acids surrounding each hit, such as one, two, three . . . 
ten or more, are replaced. 

The nucleic acid molecules of SEQ ID Nos. 563-725 and the 
encoded proteins (SEQ ID Nos. 1-562 and 726-728) are also provided. 
Recombinant AAV and cells containing the encoding nucleic acids are 
provided, as are the AAV produced upon replication of the AAV in the 
cells. 

Methods of in vivo or in vitro production of AAV or rAAV using 
any of the above nucleic acid molecules or cells for intracellular 
expression of rep proteins or the rep gene mutants are provided. In vitro 
production is effected using cell free systems, expression or replication 
and/ or virus assembly. In vivo production is effected in mammalian cells 
that also contain any requisite cis acting elements required for packaging. 

Also provided are nucleic acid molecules and rAAV (any serotype) 
in which position 630 (or the corresponding position in another serotype; 
see Figs. 5 and the table above). Changes at this position and the region 
around it lead to changes in the activity or in the quantities of the Rep or 
Cap proteins and/or the amount of AAV or rAAV produced in cells 
transduced with AAV encoding such mutants. Such mutations include 
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tgc to gcg change (SEQ ID No. 721). Mutations at any position 
surrounding the codon position 630 that increase or decrease the Rep or 
Cap proteins quantities or activities are also provided. Methods using the 
rAAV (any serotype) that contain nucleic acid molecules with a mutation 
5 at position 630 or within 1, 2, 3 ... .10 or more bases thereof for the 
intracellular expression rep proteins or the rep gene mutants covered by 
claims 10 to 13, for the production of AAV or rAAV (either in vitro, in 
vivo or ex vivo) are provided. In vitro methods include cell free systems, 
expression or replication and/or virus assembly. 
0 Also provided are rAAV (and other serotypes with corresponding 

changes) and nucleic acid molecules encoding an amino acid replacement 
by N at Hit position 350 of AAV- 1, AAV-3, AAV-3B, AAV-4 and AAV-6 
or at Hit position 346 of AAV-5; by I at Hit position 462 of AAV-1, 
AAV-3, AAV-3B, AAV-4 and AAV-6 or at Hit position 458 of AAV-5; by 
5 either R, L or Y at Hit position 497 of AAV-1, AAV-3, AAV-3B, AAV-4 
and AAV-6 or at Hit position 493 of AAV-5; by N at Hit position 517 of 
AAV-1, AAV-3, AAV-3B, AAV-4 and AAV-6 or at Hit position 535 of 
AAV-5; by S at hit position 543 of AAV-1 and AAV-6 or at hit position 
542 of AAV-3, AAV-3B and AAV-4 or at hit position 561 of AAV-5; by S 
at hit position 549 of AAV-1 and AAV-6 or at hit position 548 of AAV-3, 
AAV-3B and AAV-4 or at hit position 567 of AAV-5; by either D or S at 
Hit position 599 of AAV-1, AAV-4 and AAV-6 or at Hit position 600 of 
AAV-3 and AAV-3B; by P at Hit position 602 of AAV-1, AAV-4 and AAV- 
6 or at hit position 603 of AAV-3 and AAV-3B or at hit position 589 of 
AAV-5 in order to increase Rep proteins activities as assessed by AAV or 
rAAV productivity. Methods using such AAV for expression of the 
encoded proteins and production of AAV are also provided. 

Since modifications will be apparent to those of skill in this art, it is 
intended that this invention be limited only by the scope of the appended 
claims. 
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